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FOREWORD 


The investigation described in this report was performed by ICS during 
the period from November 1980 to August 1981, under Contract Number NASl- 
16303 for the National Aeronautics and Space Administration, Langley 
Research Center, Hampton, Virginia. This work was sponsored by the Avionics 
Technology Research Branch of the Flight Electronics Division. Dr. David 
R. Downing and Mr. Wayne H. Bryant served as Technical Representatives 
monitoring this contract.. 

The Princeton University flight crew who participated in the NAVION 
aircraft flight tests were Barry Nixon, the safety pilot, and George K. 
Miller. Wayne Bryant was the evaluation pilot for Langley Research Center. 
The author was aided in PIF flight control software development by Joseph 
Old, with the Research Triangle Institute, and Wayne Bryant. The author 
wishes to express his appreciation to Dr. Nesim Halyo of ICS for many 
stimulating discussions about aircraft control and autopilot design. 

Use of tradenames or names of manufacturers in this report does not 
constitute an official endorsement of such products or manufacturers, 
either expressed or implied, by the National Aeronautics and Space Admini- 
stration. 
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INTRODUCTION 


Commercially available General Aviation (GA) autopilots are currently 
in transition from an analog circuit system to a state-of-the-art computer 
implemented digital flight control system (Ref. 1) . Advantages of the 
digital autopilot are well known. These advantages include enhanced modes, 
self-test capacity, fault detection and greater computational capacity. A 
digital autopilot's computational capacity can be used to full advantage 
by increasing the sophistication of the digital autopilot's chief function, 
stability and control. Direct digital design techniques along with proven 
computerized design and evaluation tools should be used to efficiently 
develop low-iteration rate advanced digital autopilots. 

The General Aviation Terminal Area Operation Research (GATOR) program 
of the National Aeronautics and Space Administration is developing and 
evaluating advanced flight control and display concepts that make use of 
the recent advances in digital flight control theory, digital control hard- 
ware and electronic displays. The GATOR program uses the Princeton Avionics 
Research Aircraft to flight tests digital controllers programmed in the 
ROLM 1666 digital flight computer. Figure 1 shows the flight tests 
operation facilities. Recent work in the GATOR program evaluating electro- 
nic displays has been reported in Ref. 2. 

The purpose of this report is to present the development, design and 
flight test results of a proportional-integral-filter (PIF) digital auto- 
pilot design flown on the research aircraft. The PIF control law is a 
direct digital design operating at a low iteration rate (10 cycles per sec) 
that employs modern multivariable control theory to compute control gains . 

A basic block diagram of the control law is shown in Fig. 2. The PIF 
control law is specially structured to be designed using linear models but 
implemented to control and stabilize non-linear aircraft dynamics. The 
PIF control law does not require trim information in flight. The PIF 
control law has features to accommodate computation delay, suppress high 
frequency noise and track commands using integral control. An earlier 
version of PIF has been investigated and successfully flight tested in the 
VAIiT autoland program (Refs. 3 and 4) . The new feature is the command 
model shown in Fig. 2 which is different for each autopilot mode. 

A. AUTOPILOT SUMMARY 

The General Aviation 3-axis PIF autopilot incorporates recent 
theoretical developments in multivariable command generator tracking to 
easily design roll select, pitch select, heading select, altitude select 
and glideslope/localizer capture and hold autopilot modes. The roll and 



pitch select autopilots hold the aircraft attitude that exists at engage 
and can track pilot input from a pitch wheel and roll turn knob. In the 
heading select mode, the pilot enters in the desired heading, the auto- 
pilot command model commands the PIF control law to bank the aircraft, fly 
a coordinated turn in the shortest direction to the new heading, and roll 
the aircraft out of the turn to smoothly capture the new heading. In the 
altitude select mode, (which currently exists in abbreviated form in only 
the most sophisticated and expansive commercial GA autopilots) , the pilot 
enters in the desired altitude and the autopilot command model commands 
the PIF control law to gradually ascend (or descend) the aircraft, establ- 
ish a constant rate of climb (or descent) then smoothly capture and hold 
the new altitude. The glideslope/localizer autopilot modes use command 
models to construct horizontal and vertical paths the aircraft flys to 
capture and track precision Instriament Landing System (ILS) beams. These 
approach modes have a built-in potential to operate under and take advant- 
age of the microwave landing system (MLS) . 

The PIF autopilot uses a standard set of General Aviation instrument 
sensors; rate gyros, accelerometers, attitude gyros, barometeric altimeter, 
ILS beam receiver and airspeed sensor. Analog prefilters and digital 
complementary filters are used to smooth sensor output for improved ride 
quality. 

B. PIFCGT COMPUTER PROGRAM 


A second purpose of this investigation is the development and docu- 
mentation of an easy-to-use computer aided design program for PIF. A PIF 
autopilot design computer program (PIFCGT) which can accommodate most GA 
aircraft dynamics is described in Ref. 5. 

All the designs described in this report are obtained with the PIFCGT 
program. The PIFCGT program uses the NASA Langley developed ORACLS sub- 
routine package. Ref. 6, to solve matrix equations encountered in modern 
control design. The PIFCGT program internally constructs linear aircraft 
models from user supplied nondimens ional stability derivatives. The user 
design parameters are the quadratic weights in the Linear Quadratic Reg- 
ulator Cost function. Evaluations are performed using closed-loop equi- 
valent s-plane eigenvalues and simulations using linear aerodynamics. 

C. DISCUSSION OF RESULTS 


All autopilot modes were successfully programmed and successfully 
flight- tested using a sampling interval of 0.1 seconds. Flight tests at 
different airspeeds were performed in heavy, medium and light turbulence. 

The PIFCGT program was used to successfully adjust and completely redesign 
some autopilot modes during flight testing. In some instances flight test- 
ing, redesign, flight code reconfiguration and flight testing of the modi- 
fications occurred daily. Optimal control quadratic weight selections 
based on in-flight pilot comments of preference among different gain sets 
were performed using this procedure. An analog prefilter study was performed 
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using different filter time constants to suppress aliasing effects. An 
analog prefilter time constant of 0.1 sec gave a good compromise between 
noise suppression and closed-loop stability. 

D. ORGANIZATION OF THE REPORT 


There are a number of equations in the report and an attempt is made to 
clarify the importance of each equation. An equation with "( )" around the 

equation number is a standard result. An equation with " [ ] " around the 

equation niimber is an important equation useful in understanding and imple- 
menting the PIF control law. An equation that is in a bordered region (such 
as the PIF equations - Eqs. 56 - 58) is an equation that was actually pro- 
grammed in the onboard computer. 

Chapter II presents 6in overview of the aircraft, PIF control law and 
flight computer. A summary of advanced commercially available autopilot 
functions is also given. Chapter III is a detailed mathematical discussion 
of the theory behind the PIF control law. A complete derivation of the control 
law is presented. The chapter ends with a block diagram of the PIF software 
organization programmed in the ROLM computer. The command model parameters, 
control gains, filter gains and design quadratic weights are detailed in 
Chapter IV . Time histories from linear simulations show the command model and 
aircraft response for different autopilot modes. Chapter V discusses the flight 
tests results. Chapter VI summarizes the report and suggest areas where further 
developments not within the scope of this work could improve the PIF autopilot. 
Research results which advance the state-of-the-art of the linear quadratic 
regulator (Refs. 7, 8 and 9, for example) advanced the state-of-the-art of the 
PIF control law. Linear and nonlinear aircraft equations of motion are pre- 
sented in Appendix A . The NAVION aircraft aerodynamic model and stability 
derivatives used in this report are listed in Appendix B . The linear and non- 
linear command generator models and command control systems are described in 
Appendix C . The relationship between feedforward control and trim is reviewed 
in Appendix D . Appendix E discusses discrete structures useful in understanding 
how the PIF discrete plant representation is constructed. 
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II. CONTROL SYSTEM OVERVIEW 


A. VEHICLE, SENSOR, AND CONTROL CHARACTERISTICS 


The PIE autopilots are designed for the Princeton Avionics Research Air- 
craft (ARA) . The Princeton ARA is a NAVION airplane and is typical of a 
complex single engine general aviation aircraft. Independent control of the^ 
aircraft is provided by commands to the elevator, ailerons, rudder, throttle, 
and direct-lift flaps. The autopilot designs are restricted to use elevator, 
ailerons and rudder. These fast acting control surfaces are driven by hydraul- 
ic servos orginally fitted to the B-58 aircraft. 

Characteristics of the control actuators are summarized in Table 1 (taken 
from Ref. 10). The actuator dynamics are not included in the PIF design 
model. The time constants in Table 1 presume the actuators can be modeled by 
a first order lag. The autopilot designs are allowed full control authority 
displacement. The maximum specific moments are given at 53.8 m/s airspeed. 

The aerodynamic model for the NAVION is given in Appendix B. The estimated 
stability derivatives in the model are only available at one flight condition 
(44.0 m/s airspeed, straight and level flight) in the aircraft's normal oper- 
ating speed range (33.3 m/s (65 kt) to 69 m/s (135 kt) ) . All the autopilot 
control laws are designed using the aerodynamic model at the one available 
flight condition. This restriction is typical of the type of aircraft model 
information available to commercial general aviation autopilot designers. Step 
response comparisons in Figs. 3 and 4 between the linear model and the aircraft 
response also shows the linear model only "approximately" matches the aero- 
dynamic response, another typical autopilot design uncertainty. 

Neglecting disturbance inputs, the nonlinear equations which govern 
vehicle motion are presented in Appendix A, Eqs. 73 to 76. The total aircraft 
states and controls can be divided into nominal and perturbation components 


X 

+ 

Ax 

[1] 

— o 



u 

+ 

Au 


— o 



[2] 


The linear time-invariant model of the perturbation states of the aircraft is 
derived from the nonlinear equations at the design flight condition as shown 
in Appendix A and has the following form 
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[31 


Ax = A Ax + B Au 

Aj£ is the rigid-body state vector of the aircraft 
T 

Ajc = [ Au Aw Aq A9 Ax Az Av Ar Ap A(j) A\p Ay] 


and Au is the control vector of the aircraft previously discussed, 
A_u*^ = [A6g^ A6g^ A6^ ] 


[5] 


The NAVION research aircraft is fully instrumented and all sensor outputs 
can be recorded at the digital control law sampling interval. The PIF auto- 
pilot designs are restricted to only use sensors commercially available in 
typical avionics packages. These sensors include angular rate gyros for all 
three axis, lateral and vertical body- axis accelerometers, attitude and heading 
gyros, indicated airspeed, barometric altimeter, localizer/glideslope angular 
deviation measurements and an indication of when the aircraft crosses the outer 
marker (see Chapter V, Section D) . Measurement noise standard deviations and 
bias standard deviations estimated from in-flight recordings are shown in 
Table 2 for these sensors. Sensors available, but not explicitly used in the 
design, include angle-of-attack and sideslip vanes, radar altimeter, and range, 
azimuth and elevation from a two-way ground-to-air digital data link to radar. 
The autopilot configuration, command model designs and feedback gains would 
change if these sensors were made available. 

The general arrangement of the ARA system is shown in Fig. 5. The air- 
craft is flown by a two-man crew during all research; a safety pilot and an 
evaluation pilot. The safety pilot can disengage the control system at any 
time. Disengage software is also included in the control software for some 
flights. The control/display panel is illustrated in Fig. 6. The interface 
between the evaluation pilot and the autopilot is made through the keyboard/ 
display shown at the bottome center of Fig. 6. 

B. ROLM COMPUTER CHARACTERISTICS 


The ROLM 1666 computer is used to implement the PIF control law. The 
ROLM 1666 is a state-of-the-art digital flight computer designed to MIL-E-5400. 
The control law is implemented in a high level language (FORTRAN) , operates in 
floating point arithmetic, and can be modified easily using text editors. 

The existence of these features permitted efficient and rapid implementation 
(1 week for ALT SEL, HDG SEL, ROLL SEL, and PITCH SEL, 1 week for APR GS and 
APR LOC) and evaluation of the autopilot modes. Table 3 presents the main 
features of the ROLM 1666 (taken from Ref. 10) . 
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C. • SAMPLE SURVEY OF COMMERICALLY AVAILABLE AUTOPILOTS 


A sample survey of advanced 3-axis commercially available autopilots 
is shown in Table 4. The results in the survey are taken from autopilot 
manuals, Refs. 11 - 16 and a NASA Langley Sponsored GA Survey report. Ref. 

17. Also shown in the last row in Table 4 are the PIF controller modes whose 
designs are discussed in this report. Most of the autopilots shown in Table 
4 have become available within the last three years. Some show trends in 
the direction of full digital implementation. 

An autopilot is an add-on feature to GA aircraft. The units in an auto- 
pilot typically include actuator servos, sensors, mode controller, displays 
and an autopilot computer. The mode controller is the interface between the 
desires of the pilot and the available features of the autopilot shown at the 
top of Table 4. Autopilots are usually not command augmentation systems. In 
most operations the pilot pushes a button to activate the desired mode, then 
the autopilot holds the current flight conditions of the aircraft ( 9 , <p , ijj , 
z etc. ) . The capability to maneuver the aircraft is available in advanced 
autopilots through a select feature. The pilot enters a new desired condition 
(altitude, heading, pitch) using alphanumerics, heading bug or wheel and the 
autopilot commands the aircraft surfaces to cause the aircraft to satisfy the 
new flight condition. 

How most commercial autopilots are designed and implemented in hardware 
and software is not public information. One available source for autopilot 
design techniques is Ref. 18. Results in Ref. 18 indicate that commercial 
autopilots tend to be nested regulator designs. A pitch hold mode, when 
activated, forms a regulation error between the pitch angle at the time of 
engage and the pitch angle measured thereafter by the attitude gyro. The 
error is fed back through a gain to the elevator. The gain is easily design- 
ed by trial and error. Other longitudinal states are rarely fed back for 
improved stability unless the aircraft is aerodynamically deficient in some 
respect ( poor short period damping, etc.). The economics dictate using few 
sensors. An altitude hold mode is nested into the pitch hold mode by forming 
an altitude command error and using the error to adjust the pitch command. 

A similiar operation is performed with IAS hold. Nested designs usually 
retain control feedback gain values for each loop even after engaging the 
next control feedback loop in the hierachy. References 11-16 indicate that 
altitude select in advanced autopilots is performed through a series of auto- 
matic mode switches: vertical speed hold to altitude capture to altitude 

hold. 


The lateral directional modes have a similiar nested design. Roll angle 
error is fed back to aileron. A yaw damper is used to feedback washed-out 
(high pass filtered) yaw rate to rudder to improve Dutch roll damping. A 
nested heading hold is performed by feeding back heading error to the roll 
command. An additional outer loop is made during landing approaches by feed- 
ing back localizer deviation error to the heading command. 

D. PIF CONTROL LAW 


The PIF control law uses a different design approach from commercially 
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available autopilots. A siaitunary of control law features are as follows: 

• The multivariable PIF control law is designed using modern control 
theory and is a direct digital design. 

• The control law is composed of five basic units: a command generator, 

state feedback for stability, command error integration for command 
tracking, complementary filtering for sensor noise suppression and a 
control filter for computation delay accommodation as shown in Fig. 2. 

• The control law is structually formed to stabilize and command aug- 
ment any nonlinear systems. PIF works particularly well for aircraft 
nonlinear dynamics as has been demonstrated by flight tests and non- 
linear simulations for differing aerodynamic vehicles (PIF-Tandem 
Rotor Helicopter, Ref. 3, Pi-High Performance Fighter, Ref. 19). 

• The control law is designed using perturbation variables but does not 
use nominals, Uq and Xo/ the implemented control logic. Computing 
nominals in flight as suggested in Ref. 20 is impractical for GA auto- 
pilots and not necessary because the control law uses the incremental 
form. Ref. 21. 

• PIF is a command augmentation system design. To perform the "select" 

feature a command model representing simplified aircraft dynamics 
uses a nonlinear control law to transfer internal command states from 
one desired condition (ex: altitude) to a different desired condition. 

The PIF control law tracks the command model states . The command 
tracking is performed by combining the derivation of the PIF control 
law with the theory of feedforward control. Refs. 22 and 23. The 
derivation is detailed in Chapter III, Section B. The command 
models are discussed in Appendix C for each mode shown in Table 4. 

• A computer program is available to design any of the autopilots in 
this report for any aircraft with known stability derivatives. The 
program automatically linearizes the nonlinear aircraft dynamics to 
construct A and B in Eq. 3, computes feedback and feedforward gains, 
determines closed- loop eigenvalues and simulates the completed PIF 
design. 

• PIF is not a nested control design in implementation. When a new 
command mode is engaged, the PIF gains are updated to optimally 
match the new command mode. Alternate software sections of code are 
executed for each command mode. Approximately 75% of the PIF control 
software, however, is always executed independent of the command mode. 
Nested control law implementation is to be distinguished from a nested 
construction of the command error. The increase in complexity of the 
command error from Eq. 142 from the ROLL SEL mode, then Eq. 145 for 
the HDG SEL mode to Eq. 191 for the APR LOCR mode is a result of 
nested construction. 
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III. DIGITAL FLIGHT CONTROL SYSTEM 

This chapter derives the PIF control law. A derivation of the PIF 
control law optimally transferring a linear plant from tracking one constant 
command input to another is given in Ref. 3. This chapter generalizes the 
derivation to include feedforward control and model output tracking and is 
unique to this report. The derivation is accomplished by formulating a 
continuous- time optimal regulator problem, converting the problem to a sampled- 
data regulator, incorporating the "star trajectory" (defined later) , solving 
the discrete-time regulator problem, then constructing the final implementable 
solution. 

The PIF drivation uses assumptions about the command model control input 
and incorporates knowledge of digital structures during the construction of 
the aircraft design model. These assumptions and structures are reviewed in 
the next section and Appendices D and E. 

A. COMMAND MODEL AND FEEDFORWARD CONTROL 


The main objective of the PIF control law derivation is to cause the out- 
put of a linear aircraft model given by 


= H Ax + D Au 


[ 6 ] 


to optimally track the output of a linearized command model. 


Ax = A Ax + B Au 
^m m — m m -m 


[7] 



H Ax 
m — m 


+ D Au 
m — m 


[ 8 ] 


The PIF derivation causes A^ to track A^ assuming the command model control 
input, Ai^, is constant for t ^ 0. The model control input is not constant 
in implementation. The rest of this section explains the reason behind the 
A^ assumption and its impact on control system performance. 

There are at least three methods, or approaches, to solving the command 
model tracking problem so that Ay tracks Ayjjj. 

1) Ignore the command model dynamics, assume Ay^ is constant and feed- 
back the command error A^ - A^j^ with a high gain feedback loop. 

2) Include a judicious part of the command model, assume A^ is constant 
and feedback the command error A;^ - A^m with moderate feedback so 
that A;^ tracks A^ when Aum is constant. 
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3) The exact command model is used, is assumed to be arbitrary, 

after a brief initial condition mismatch transient, A^, tracks A^ 
perfectly no matter how Au^^ changes. 

The first approach has been used successfully in many practical designs. The 
PIF control law used in the VALT program. Ref. 3, is designed using the first 
approach. Velocity commands are tracked by the PIF control law and are assumed 
to be constant during the design. In implementation the velocity commands have 
low frequency variations caused by feeding back position error as the helicopter 
is required to track a complex 4-dimensional trajectory. 

A state-space solution to the third approach is given in Ref. 23 and 
is similiar in formulation to exact model matching (Refs. 24 and 25). Ref- 
erence 26 shows that certain adaptive controllers are adapting their internal 
gains to achieve exact model matching in steady state. Experience with Approach 
3 has shown that for autopilot designs, the approach is impractical. Exact 
pole/zero cancellation is required, the control law is sensitive to noise, the 
required effector response is not achievable by motor driven rate limited actu- 
ators and the discrete-time solution requires future values of Aujjj. 

A practical method and the one pursued in this report is method 2. When 
Ax^ is changing, and became mismatched and their error is governed by 

the control law in a manner similiar to approach 1. When A uj^ is constant, A^ 
eventually tracks A^ as in approach 3 but suffers none of approach 3's dis- 
advantages that stem from the arbitrary Ax^ assumption. The PIF control law 
using approach 2 is designed to optimally govern the discrete time transient 
response of A^ when Ax^ changes once and remains constant thereafter. In 
practice the control law performs well if Ax^ changes intermittently or 
slowly varies. 

The feedforward control solution where Ay tracks Ai^ with A_^ constant is 
presented in Appendix D. The feedforward control solution is independent of 
the control law used and is as follows: 


1 



Ai2 


Ax , 
— m,k 

1 
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CM 

^22 


Au 
— m 


[9] 


The output A^ tracks A^ when Ax follows Ax* and the feedback control is A^* , 
i.e., the state and control star trajectory. The solution is unique if the 
nxmber of plant outputs equals the number of plant controls. The feedforward 
matrices A^j satisfy Eq. 227; a solvable matrix algebraic equation (Ref. 5) . 

The star trajectory is a convenient notational abstract and is not 
physically generated when the control law is implemented. All the linear 
simulations of the PIF control law (see Chapter IV, Section D) however, show 
the star trajectory for clarity. The relationship between the linear star 
trajectory and trim in Eqs. 228 and 229 is assximed to remain valid and is 
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exploited when the PIF control law is derived in the next section. 
B. PIF DERIVATION 


The continuous- time design model and cost function for PIF are as follows 
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( 10 ) 


( 11 ) 


The linear model of the aircraft shown in Eq. 3 is augmented with an equation 
for control rate and an equation for the integration of the command output, 

A^. The augmented plant model and quadratic cost function are in standard 
infinite-time linear quadratic regulator form. 

Four problems a digital flight control system should address are (1) the 
type of control data-hold employed, (.2) the effect of the computation delay, 

(3) high frequency noise suppression and (4) trim accommodation. These 
problems are discussed in Appendix E. The PIF control law addresses the first 
three problems by including Au dynamics in the model and weighting Av in the 
cost function. The trim problem is partially addressed by using the inte- 
grator, A^, in the augmented model. 

Precisely converting Eqs. 10 and 11 using the sampled-data regulator. 

Ref. 27, where Av^^ is assiamed constant between equally spaced sampling points, 
yields a preliminary version of the discrete-time PIF design model 
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The matrices (j), T, and are defined in Eqs. 218, 219, 239. The discrete- 
time model in Eq. 12, if used to design a PIF control law, requires that the 
control law be implemented with a triangular data-hold (compare Eqs. 12 and 
238) . Similarly, the discrete dynamic model for the integrator in Eq. 12 is 
a discrete representation of the implementation 
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> t i t 


(13) 


A£(t) = A^(t) dt 

k 


k 


The required analog integration of Ay(t) between samples as shown in Eq. 13 
is not desired in implementation. 

The problems with using Eq. 12 as the discrete model representation are 
eliminated using approximations. The discrete cost function obtained using 
the sampled-data regulator and shown later in Eq. 22 is not changed when the 
approximations are introduced. Hence, the optimization problem PIF solves is 
not precisely a sampled-data regulator but the effect is small, (Ref. 28) , 
because the terms eliminated are of the order At^. 

The first approximation is performed by setting to zero rn Eq. 12. 
When is zero, the required implementation for Ai^ changes from a triang- 
ular data hold to a zero-order hold. The resulting control law is also re- 
quired to be implemented with a modeled one sample period computation delay 
(compare Eq. 12 with F^ zero and Eq. 243). 

The simplification of the integrator in Eq. 13 can be performed by using 
any number of discrete integration formulas, such as 



= A?^ + At A^j^ 

Euler 

(14) 


= + x'"2k+l "^k> 

Trapezoidal 

(15) 

"4+1 

' "4 "A^k+i 

Rectangular 

(16) 


Euler integration is used for simplicity. Euler integration implies that the 
integration error is delayed two sampling intervals in implementation (see Eq. 
51) . Using rectangular integration instead of Euler integration would cause 
the integration error to be delayed only one sampling interval but would pro- 
duce difficult but not insurmountable implementation problems. 

The star trajectory in Eq. 9 is introduced into the PIF model by defining 
the variables 


Ax, 
— k 


Ax, - Ax,* 

— k — k 

Au^ - 

Av^ - iv* ^ 


(17) 

(18) 
(19) 
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A4 = Aik - * 


( 20 ) 


The star trajectory for is defined later (see Eg. 37) . Combining the 
simplifications with Eqs. 17 to 20 results in the discrete PIF model,. 
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[ 21 ] 


The integrator model in Eq. 21 is obtained using the following steps: (1) A^ 

in Eq. 14 is replaced with A^ - Ay ^ ^ ^ , (2) ^ is replaced with hAxj, + dA^ 
using Eg. 6, (3) A^^]^ is replaced with Ay]^ using results from Appendix D, 

(4) A^ is replaced with hAj^ + dA\^, (5) A^* is subtracted from both sides of 
the integrator equation to obtain the model for A|j^, (6) the model for A|j^ 
replaces the model for A^ in Eg. 12 to obtain the model following integrator 
dynamics. 


The cost function is the sampled-data cost function obtained using Eqs. 
10 and 11 but starting at k = -1 instead of k = 0. 
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The cost function starts at k = -1 since Au^ and Av^_j_ are unspecified. The 
cost function is minimized by assuming changes at k = 0 and remains 

constant thereafter. The plant is assumed to be tracking the star trajectory 
for = Ai ^^-2 = ••• previous to the time instant k = 0. 

The changing model input at k = 0 introduces a contradiction in the 
value of Axq. The ideal plant dynamics indicate 


Ax* = $Ax*^ + FAu*^ 
where 
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( 24 ) 
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At k = 0, a different values for Ax* is obtained using Eq. 9 which uses 
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Substituting Eq. 24 into Eq. 23 and simplifying does not result in Eq. 25. 

The contradiction can be resolved in the model and feedforward equations , 
by advancing the time index on Au one step as follows 
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(26) 

(27) 


Siibstituting Eq. 27 at k = -1 into Eq. 23 and simplifying produces 

Axi = (j>Ax* + TAu* = ^Ax +A, „Au „] + T[A„^Ax ^+A„„Au _] 

—0 — -1 — _1 11 — m,-l 12 — m,0 21 ^,-1 22 — nijO-* 


A, T Ax ^ + A. Au ^ 
11 ^,0 12 -m,0 


(28) 


the required result. The derivation of Eq. 28 uses Eq. 227 to perform the 
required substitutions. The redefined model does not alter the derivation 
in Appendix D for k > 0, but does beneficially change the PIF implementation 
and command response. 

Equations 21 and 22 are in the proper regulator format and yield the 
solution, 


(Au^ - ^^k-1^ " “^^k-l ” ^1 ^\-l ^2 ^-k-1 ^3 ^^k-1 (29) 

where the gains C^, and determined from the algebraic Riccati 

equation solution'^associated with the optimal regulator (Ref. 5) . The inte' 
grator star trajectory is still unknown in Eq. 29. 

It is well known that the minimum cost for the optimal regulator is 
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min J 


( 30 ) 


' >’*20 


where P is the steady-state solution to the discrete algebraic Riccati 
equation. It follows that the minimum value of the perfoianance index, shown 
in Eq. 22, is (for fixed Aj^*) 


min J 


P P 





(31) 


Further minimizing J by setting 9J/8A^_-L = 0, because A^* is unspecified 
produces 





(Au - Au* ) + p'^^CAx 
-1 ~-l --1 


Ax* )J 


(32) 


Recalling that Au_j^ and Ax_-j_ are assumed to be tracking the star trajectory 
for A^ produces 


Ax - Ax* = (A Ax , + A, „Au -) - (A Ax + A, „Au ) 
—_l 11 — m,-l 12 — m,-l 11 -m,-l 12 m,0 


(33) 
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S\ibstituting Eqs. 33 and 34 into Eq- 32 and simplifying produces 

(35) 

(36) 
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is assumed to be tracking a star trajectory for A^,_i previous to k = 0, 
hence, Eq. 35 implies that the optimal star trajectory for the integrator state 

is 


Ag* = [0 A] 


Ax , 

— in,k 

Au , 

— tn,k+l. 


k — i» • • • 


The control law in Eq. 29 can be written as follows 

Au, == Au, , + At Av, 


-k 


-k-1 


-k-1 


Av^_l = + ^2 ^-k-1 ^3 ^^k-1 


^4 = A4_3_ + At A4_^ 


(37) 


(38) 

(.39) 

(40) 


Subtracting Avj ^_2 from Avi ^ _ t in Eq. 39 changes Eq- 39 to 


^4-1 = (I + At C^) Av^_ 2 + C^(Axj^_^ - Ax^_2> 

The perturbation variables are eliminated from Eqs. 38 and 41 using a large 
number of substitutions and cancellations. An example of the substitutions and 
cancellations is demonstrated for the incremental state feedback term 

C2(A4-i “ A4-2^ 

The state feedback increment is simplified as follows 


^\-l ■" '^■^-2 ^^-^-1 " ^-^-2^ " ^'^-k-1 ^-*k-2^ 

^-1 “ ^-2 " ^-k-1 ■ -k-2^ “ ^^,k-l ' -0,k-2^ 

^^k-1 - ^^k-1 = 4l^^^,k-l " ^^,k-2> 42^^^, k ■ ^^,k-l^ 

^^,k-l " ^^,k-2 " ^^,k-l " ^,k-2^ ■ ^^,0,k-l " ^,0,k-2^ 
^^,k ^^,k-l ^^,k ^,k-l^ ^^,0 -^,0^ 


(42) 

(43) 

(44) 

(45) 

(46) 
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but 


1 1 “ X- , _ - Ait(x _ , 1 - X „ - „) + A-„(u ^ - u ( 47 ) 

0,k-l 0,k-2 11 -^,0»k-l “tn,0,k-2 12 — m,0 — m,0 


hence 


Ax, - - Ax, ^ = (X, T - X, _) - A-, (X , 1 - X , ry) 

— k-1 — k-2 — k-1 — k-2 11 -m,k-l -^,k-2 


- A, „ (u , - u , , ) 

12 -THjk — m,k-l 


(48) 


Using logic similiar to Eqs, 42 through 48, the PIF control equations reduce 
to the following 


u, = u, , + At V, T + A„, (x , - X , ) 

^ -k-1 —k-1 21 — m,k — m,k-l 


e, 1 = X, , - A, , X , , 
— k-1 — k-1 11 — m,k-l 


Vfc.i = (I + At Cj) v^_2 + ■ -£k-2^ S<2k-2 


+ ( Cj A 


(49) 

(50) 


^,k-2^ 


(.51) 


In all of the autopilot models discussed in this report, ~^2\ ^ zero 

matrix. When h .21 is eliminated from Eqs. 49 and 51, the only changes to the 
PIF control law over Ref. 3 are the propagation of the command model to generate 
^ new matrices are used to compute the feedforward matrix for (i^ 
in Eq. 51 and (^_i ~ .Sk-2^ used for state feedback instead of ( 2 ^-x “ Xfc-2^ ■ 

The last modification needed for the PIF control law is due to the sensors. 
The mathematical model of the aircraft in Eq. 10 uses the body axis perturbation 
states Au, Aw, Aq, A0, Az, Av, Ar, Ap, A(j>, Aifj, Ay, which are related to the 
perturbation sensors output Az^ = [Av, Aa^/ Aq, A0, Az, Aa^., Ar, Ap, A(p, Ail>, Ay] 
through the relationship. 



[52] 


The matrices Cj^ and are obtained using Eq. 135 for Aa^ and Aay, and Eq. 137 
for Av. The rest of the elements in are ones in the appropriate place along 
the diagonal for the body axis states that are directly measured. Equation 52 
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is used to incorporate sensor feedback in Eq. 39 by using the following 
substitution 


'1 '=ul‘ 

UiJ Lo iJ La^J Lau^J 


The feedforward matrices are similarily altered 


fc cirAx*"! fc cir 
X I I = ^ ^ 

_o I. Lau*J Lo iJ ^ 


A- T A. „ r Ax 
11 12 — m,k 

A A J _Au , 

21 22 -^,k.+ 


-^91 ^79- -^^.V+1- 


J L^21 ^22-' 

Eliminating A 21 and using sensor feedback, the final implementable form for the 
PIF controller advanced one time step is 
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[61] 


The PIF simulation program in Ref. 5 and the flight computer software discuss- 
ed in the next section use Eqs. 56, 57 and 58 for implementation. Figure 7 
shows a detailed block diagram of the PIF equations which includes the command 
model equations. 

C. FLIGHT COMPUTER SOFTWARE 


The filter and PIF control logic implemented on board the ROIiM 1666 flight 
computer is presented in this section. The flight computer software uses a main 
executive routine which calls subroutines for each part of the control computa- 
tions. The flight software is written strictly to test the autopilot modes. 
Safety features, sensor failure detection, actuator failure detection, bad data 
editing (except for radar and yaw gyro) , graceful degradation features and 
built-in-test are not included in the software. 

The sequence for a single control computation frame is illustrated in the 
fianctional flow diagrams of Figs 8 to 11. The operations performed by the 
blocks in the functional flow diagram are described next: 

MAIN EXECUTIVE 

ONE TIME TASKS - Communication channels are open, and logical variables 
are initialized. 

TERMFLEX - Software for a push button panel is executed. The alpha- 
numeric display prompts the pilot concerning the desired command mode 
and command inputs. The pilot inputs are stored and logical variables 
are initialized. 

READ A/D SENSORS - Sensor output from the analog to digital converters 
(which is synchronous with the 10 per second computation loop) are scaled 
and stored in common blocks. 

FILTER INITIALIZATION - The filter variables are reset to the sensor 
measurements or zero until the pilot engages the filters (depressing 
INITIAL) . 

EXECUTE FILTERS - The barometric complementary filter and yaw complementary 
filter computations are executed. 

CONTROL INITIALIZATION - Control increments and command errors are set to 
zero and control commands are updated to actuator positions until the pilot 
engages the control law (depressing ENGAGE) . 

CONTROL - The PIF control law is engaged. The command values are the 
flight conditions that exist at the time of engagement. New commands 
obtained in the TERMFLEX subroutine are activated by the pilot in flight 
(depressing COMMAND) . 
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UPDATE SERVOS - The control actuator commands are prepared to be released 
at the beginning of the next cycle. 

WAIT REST OF INTERVAL - Execution pauses until the beginning of the next 
cycle. 

CONTROL 

COMPUTE - The control law computations are performed unless this is the 
first pass or the pilot has pressed the COMMAND button. 

NEW COMMAND - When the pilot presses the COMMAND button the desired com- 
mand change is computed and command model following begins execution, 

FIRST PASS - A two-dimensional array contains the PIF control gains 
stored columnwise for each autopilot mode. In the first pass the desired 
gain set from the two-dimensional array is loaded into a vector for later 
use in FPIFl. 

APRFILT - Estimates of lateral position, perpendicular distance to the 
glideslope and range are computed in this subroutine. 

APR GS, ALT SEL, PITCH SEL, APR LOG, HDG SEL, AND ROLL SEL - These are 
the subroutines which accept command inputs, compute command model controls 
and propagate command model states to generate desired trajectories. The 
equations used in these subroutines are given in Appendix C. 

FPIF2 - The feedback and feedforward error quantities for the PIF control 
law are computed in this subroutine which is detailed in Fig. 10. 

FPIFl - The control gains are multiplied by the feedback and feedforward 
errors and the new control actuator commands are determined in this sub- 
routine which is detailed in Fig. 11. 

LOC/GS - Logic is used to test when to switch from the ALT SEL and HDG 
SEL modes to APR GS and APR LOG modes. 
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IV. DIGITAL FLIGHT CONTROL SYSTEM DESIGN 


The PIP autopilot design procedure is a sequence where the designer 
inter atively improves system performance at each step. The sequence in the 
design procedure and the parameters adjusted are shown in Fig. 12. The PIFCGT 
computer program aids the designer to perform the functions enclosed in the 
dotted block in Fig. 12. This chapter discusses each part of the design pro- 
cess and presents the design parameters for the PIF autopilot control laws. 

A. COMPLEMENTARY FILTERS AND PREFILTERS 


This section presents the designs of the complementary filter and pre- 
filters used to smooth sensor output. Originally the PIF control law for the 
NAVION was designed and flight tested with no filtering of the digital sensor 
data. Flight tests indicated filters were needed for the barometric altimeter, 
yaw gyro and lateral accelerometer. 

The barometric altimeter signal is very noisy due to quantization, atmos- 
pheric effects, and sensor performance- An example of the noisy signal is 
shown at the top of Fig. 13. A complementary filter discussed and designed in 
Ref. 3 is used to smooth the vertical height measurement. The complementary 
filter uses steady-state Kalman filter theory to determine the filter gains 
based on sensor noise statistics reported in Table 2. The 3rd order comple- 
mentary filter consists of the update equation 
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The states in the filter are vertical height (positive up) , z, vertical velo- 
city, z, and accelerometer bias in body axis b. The accelerometer output, a^/ 
in flight is biased to read zero (positive up) in straight and level flight 
and must be corrected using gravity, g. 


a ,= a ,+g* (1.0 - cos 9, cos cj), ) [64] 

z j z y k iC cC 


h]g is the sampled output of the barometric altimeter. In flight performance 
of the filter presented in Fig. 13 shows that the height estimate is favorably 
comparable to height estimated using radar. 

The yaw gyro output is from 0 to +7T and 0 to -IT. A complementary filter 
is used to smooth the sampled yaw gyro measurement when the signal crosses 
between +TT and -TT. The complementary filter design is a modification of a 
design presented in Ref. 3. First, the filter residual error is computed and 
checked for the ±7T rad crossing. 
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The propagate equation is 
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The filter gain is adjusted to decrease near ±7T rad as follows 



An estimate of using body - reference rates q and r is used in Eq. 67. The 
performance of the filter is shown in Fig. 14 where the aircraft heading crosses 
180 deg and the flight is performed without difficulty. 


All sensors began flight testing using one pole low-pass (RC) analog pre- 
filters with 0.01 sec time constants to reduce aliasing of very high frequency 
electronic noise. Initial flights with lateral accelerometer feedback yielded 
the high oscillatory response shown in Fig. 15. 7Ui experimental test replacing 
ay with sideslip feedback eliminated the oscillation. The conclusion is that 
the lateral accelerometer is sensing a moderate frequency noise, probably from 
engine vibration, which is being aliased and caused the low-frequency oscilla- 
tion in Fig. 15. The capacitor in the RC prefilter network was changed to 
produce prefilters with time constants of 0.10 sec and 0.22 sec. The flight 
test results are shown in Figs. 16 and 17. The safety pilot introduced distur- 
bances during the flight tests by pulsing the rudder pedal through the cable 
system. As the prefilter time constant increased, Dutch roll damping decreased, 
as is evident in the ay output. The 0.1 sec time constant prefilter is chosen 
as an excellent compromise between stability and noise suppression. The 0.1 
time constant prefilter was used on the ay and accelerometer outputs for 
all other flight tests in this report. 

B. AUTOPILOT COMMAND MODELS 

Each autopilot mode is determined by the desired output vector, that 

the aircraft is required to track and the model dynamics needed to change yjj, 
depending on the desires of the pilot. A complete discussion of the autopilot 
models is given in Appendix C. The model dynamics are sometimes non-linear 
and a linearized representation is needed to detemine the feedforward matrices. 
Table 5 shows the model states and aircraft states along with the controls used 
in the PIF autopilot mode designs. 

The parameter, aj., in the lateral autopilot is the roll-to-rudder gain 
used to produce a coordinated turn. The parameter varies in flight according 
to 
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a, = -0.001093 V, + 0.08268 
k k 


where is limited when 

< 39.63 THEN = 0.03937 

and > 60.93 THEN a^ = 0.01605 

The parameter a^ is determined using BETA HOLD system as discussed in Appendix C. 

The linear command models in Table 5 are simpler than the nonlinear com- 
mand models. Experience with the models has determined that including all model 
states and controls in the linear models can degrade tracking performance. The 
basic question as discussed in Chapter III, Section A, is the amount of 
derivative information needed by the design. More derivative information causes 
a faster response. If is constantly and quickly changing, as does in the 
HDG SEL autopilot, the PIF design over compensates causing overshoots, and poor 
ride quality transients, but quickly drives Y. “ "to zero. Using as u^^^ 
instead of in the feedforward control design causes the aircraft roll angle 
to lag behind the roll command. The resulting design has good ride quality 
transient response, however, as the flight tests demonstrate 


The command models use a control system to steer command model states using 
the command model controls. Each control system is presented in Appendix C. 

Each control system is designed to produce desirable paths for the PIF control 
law to track. The command model control system design parameters are shown in 
Table 6. Pilot comments are used to fine tune model control system parameters. 
The command model control system designs are independent of the aircraft model, 
hence, a PIF autopilot for a different aircraft can use the same parameters in 
Table 6. 


C. PIF DESIGN AND EVALUATION 

The PIF design involved choosing quadratic weights, evaluating the result- 
ing design using gain values, closed-loop s-domain eigenvalues and command re- 
sponse, then adjusting the quadratic weights as necessary. The PIF design pro- 
cess benefitted from the fact that during a flight test, the pilot was able to 
switch between gain sets then indicate the preferred design. The HDG SEL and 
ALT SEL autopilots were the first two designs accomplished with the pilot opinion 
approach. Once these two designs were completed, the other autopilot quadratic 
weights were easily obtained using HDG SEL and ALT SEL as a baseline design. 

The quadratic weights chosen for all the. designs are given in Table 7. 

The ROLL SEL autopilot was the first designed, did not use pilot opinions 
to obtain the weights, and was one of the last autopilots to be flight tested. 
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The ROLL SEL autopilot performs a simple function and the initial 'design did 
not have to be adjusted after flight testing. 

The other designs all have equal, continuous-time, cost function weighting 
on control rate and control position. Reference 7 shows that equal weighting 
along the control weighting matrix diagonal produces LQR designs with good 
robustness properties. 

The HDG SEL autopilot quadratic weights were adjusted, using criteria 
discussed later, until reasonable time histories and closed- loop eigenvalues 
were obtained. Different values for the quadratic weights were chosen to 
yield a baseline design and alternative designs. The pilots were satisfied 
with the stability and performance of the baseline design for the HDG SEL 
autopilot. The f light-tested control law perfomnance and linear simulation 
performance were in good agreement. 

The longitudinal ALT SEL autopilot flight tests and the linear simulation 
had only fair agreement. The PIF dontrol design program (PIFCGT) performs a 
full state feedback design without knowledge of any state filtering. The 
implemented PIF control law uses a complementary filter to estimate vertical 
height from a noisy accelerometer and barometric altimeter. 

The problem with designing the ALT SEL autopilot is shown in Fig. 18. 
Throughout the flight, atmospheric disturbances cause the autopilot to 
continually interplay height and pitch attitude in an oscillatory (phugoid 
mode frequency) manner which we classify as altitude hunting. Increasing 
the quadratic cost function weight on Az and A^^ increases the frequency and 
causes the pilot to notice the changing pitch attitude. Decreasing the weight 
on Az and A^^ eventually causes the pilot to notice the changes in height. A 
compromise is reached by increasing the weight on A^^ until the pilot does 
not notice the changes in height, then increasing the weight on A0 to decrease . 
the pitch attitude excursions. Most ALT HOLD autopilots have a "tight" pitch 
inner loop to increase phugoid damping (Refs. 18 and 29) . The effect of in- 
creasing the weight on A0 is shown in Fig. 19. When the gain set is switched 
from the low A0 weight to the design weight gain set, the pilots perceive that 
the altitude hunting is eliminated and the plane is flying straight and level. 

A difference between the two gain sets is not evident from the figures. 

The s-domain closed- loop eigenvalues for the autopilot designs are shown 
in Table 8. The s-domain eigenvalues are obtained by computing the eigen- 
values of the natural log of the discrete closed-loop system matrix as describ- 
ed in Ref. 3 for the PIF control law. 

The open- loop short period mode and roll mode are well damped and stable. 

No quadratic weights are needed on Aq, Aw and Ap. As noted in Ref. 29, another 
major problem with altitude hold is speed control. The speed mode in Table 8 
is stable, but has a large time constant determined primarily by the aero- 
dynamics, (a zero in the Az/A6^ transfer function) . Increasing the weight on 
Au greatly increases control gains and degrades altitude hold performance since 
the control law begins to use elevator to stabilize speed by diving and climbing, 
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changing altitude with the ALT SEL autopilot requires that the pilot adjust 
throttle, particularly in a climb, to maintain airspeed. The PITCH SEL auto- 
pilot speed mode is more stable, as shown in Table 8 , than the ALT SEL auto- 
pilot and uses lower gains. Closing the longitudinal PIF control design loop 
with the Az complementary filter in the loop porbably alters the closed-loop 
eigenvalues shown in Table 8 to produce the lower damped oscillatory mode 
which is evident in flight tests. 

The APR GS autopilot aircraft model with Ad dynamics is identical to the 
model with Az dynamics for straight and level flight. The difference between 
the two designs is that APR GS uses an even higher weight on A0 than the ALT 
SEL design gain set. 

Table 8 shows that the Dutch Roll damping for the lateral-directional 
autopilots is low but acceptable. A large niomber of quadratic weight itera- 
tions have been made in an attempt to reduce the ay to vgj. gain which contri- 
butes to lowering the Dutch Roll damping. No weight is placed on Av and large 
weights are placed on the angular positions. The large weighting on Ar in 
Table 7 is used to offset the zero weighting on Av. 

The reason for keeping the ay to vg^ gain low is because of the sensiti- 
vity of the PIF Gain Cg in Eq. 59 to The linear model for the lateral 

accelerometer discussed in Appendix A simplifies to 

Aa = Y Av + Y. A6 (71) 

y v 0 r 


in straight and level flight. A nonzero values for causes the matrix in 
Eq. 52 to be invertible. The parameter is absorbed into the Cg gain in 
Eq. 59. The low-pass filtering effect in PIF changes as ¥ 5 ^. varies with 
flight condition. PIF designs computed for a large number of different flight 
conditions in Ref. 3 indicate that Cg (without using accelerometer feedback) 
is essentially constant and relatively insensitive to changing flight conditions. 

A 400% mismatch between Ygj. in Ref. 30 and Yg^ identified from flight data 
using the ay, 3 , measurements and Eq. 71 was the cause for early flight tests 
of PIF to be unstable. After Yg^ was reasonably identified, no further dif- 
ficulties were encountered in flight tests due to accelerometer feedback. 

The closed- loop properties of the three APR LOG auotpilot modes are shown 
in Table 8 . The closed- loop properties of the three APR LOG autopilot modes 
are shown in Table 8 . The closed-loop eigenvalues of the three modes are 
similiar, indicating the difference between the three modes is primarily due to 
different command following properties. The quadratic weights in Table 7 are 
chosen so that the three modes have similiar closed- loop eigenvalues. 

The PIF autopilot feedback and feedforward gains are shown in Tables 9 
and 12. The gains are constant in flight and are designed at 44.0 m/s (85 kts) 
straight and level flight. The diagonal elements in the PIF gain Gg in Eq. 59, 
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based on previous studies (Refs. 3 and 19) , appears to yield a good design for 
values between 0.6 and 0.7. The control rate weights in Table 7 are adjusted 
until the diagonal elements in Cg approach this desired range. The v,5j. - v<5^ 
rudder gains in Tables 10 and 12 include the effect of Y(5j, and are actually 
in the range of 0.66 to 0.69. 

The ALT SEL and HDG SEL autopilot gains, which were obtained by pilot 
opinion between different gain sets, serve as a baseline for the APR LOG modes. 
The weightings on Ay and Table 7 were obtained by iterative adjustment 

until the y - v^^ and Y " gains were near the z - v^g gain and gave ade- 
quate tracking performance. Gain adjustment through quadratic weight selection 
for the APR LOG modes were not performed during flight testing. 

The matrix elements needed to implement Eq. 57 are shown in Table 13. 
The order of the aircraft states, x, and model states, 5^, are shown in Table 
5. The APR LOCR is the only design that requires an Sn element that is not a 
0.0 and 1.0. The 0.31 element is caused by the way (y - is constructed in 

the APR LOGR autopilot and has little effect on (y_ - since ijJ is in radians 

and y is in meters. The 0.31 gain depends on the values of kijj and ky in Eq. 191 
and the aerodynamics. 

D. PIF LINEAR SIMULATION 


The linear simulation time histories used to evaluate the PIF autopilot 
designs are shown in this section. The simulations use the linear time in- 
variant model shown in Eg. 3 to generate the time histories for Ax. The 
actual nonlinear command models are used in the linear simulation to generate 
the time histories for The ALT SEL, HDG SEL, and ROLL SEL simulations 

start with all initial conditions and commands zero. One second into the simu- 
lation a step command (altitude, heading or roll angle) is issued to the com- 
mand model system initiating control action. The APR GS, APR LOGI, APR LOCR 
and APR LOCP autopilot linear simulations begin with nonzero initial condition 
offsets in d or y position. The contorl law drives these offsets to zero. 

Each simulation plot contains two curves. The solid line curve is the 
output of the linear aircraft system model. The dotted line curve is the star 
trajectory generated by Eq. 9. The difference between the solid curve and 
dotted curve is the regulation error that is weighted in the quadratic cost 
function and forced to zero. 

The ALT SEL simulation is shown in Fig. 20 for a 30.48m (100 ft) descent 
in height. The command model deceleration is initally limited resulting in 
gradual transition to a steady descent. As the command model approaches the 
desired height, the vehicle vertical velocity is vigorously returned to zero 
resulting in a smooth capture with virtually no overshoot. The pilot prefers 
the higher acceleration at the end of the altitude change since this serves as 
an anticipatory cue that the autopilot is completing the maneuver. The height 
regulation error, which is the primary error of interest, is small; while the 
other regulation errors in Fig. 20 are large. The large regulation errors are 
caused by the slow response of the speed mode when throttle is not activated. 
Attempting to reduce the regulation errors in states other than Az degrades 
the Az response. 
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The HDG SEL mode response for a 45 degree change in heading is shown in 
Fig. 21. When the heading change is entered, the autopilot banks the aircraft 
using a constant roll rate up to a specified roll angle limit in the command 
model. The bank angle direction flys the shortest turn to the new heading. 

The command system begins rolling the aircraft out of the turn before reach- 
ing the new heading producing no heading angle overshoot. The ROLL SEL auto- 
pilot simulation in Fig. 22 executes a step command in roll angle with modest 
overshoot and rapid settling. 

The APR LOG linear simulations shown in Figs. 23 and 30 are, as expected, 
very similiar. The APR LOCI has the best capture performance. The APR LOCR 
simulation has larger yaw angle and lateral acceleration deviations while the 
APR LOOP autopilot uses larger the necessary roll angles. The simulations in 
Figs. 23 and 25 use an initial intercept angle of 0.0 degree. The flight tests 
use steeper intercept angles in the range of 30 to 60 degrees. 

The APR GS linear simulation initially positions the aircraft below the 
glideslip as shown in Fig. 26. The d position below the glideslope approaches 
0.0 as the aircraft flys towards the runway. Before the d position reaches 
0.0, the autopilot pitches down and executes a smooth capture with no overshoot. 
The distance, d, is the perpendicular distance between the aircraft eg and the 
glideslope and is discussed further in Appendix B. 

The approach autopilots exhibit highly nonlinear behavior because of the 
complexity of the command model and perform differently depending on the 
intercept initial conditions. Continued discussion of the approach autopilots 
design and flight test performance is given in Chapter V, Section D. 
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V. DIGITAL FLIGHT CONTROL SYSTEM FLIGHT TESTS 


Flight test results for the PIF autopilot designs described in the 
previous chapter are presented in this chapter. The flight tests are made 
with the aircraft and equipment described in Chapter 2. The data are from 
a series of tests that occurred over a five month period. Turbulence and wind 
conditions encountered varied from heavy to light and their effect is evident 
in the data. The plotted data includes measurements available onboard but 
not used by the autopilot (namely, radar positions, sideslip, range). When 
the autopilot is engaged, all the recorded sensor outputs show a small amount 
of bias shift as evidenced by small early transients in plots that include the 
system engage. 

The orientation of the axis system for the flight tests and for the linear 
simulation is different as discussed in Appendix B. In the flight tests, the 
vertical axis is positive up. In the linear simulation, the vertical axis 
is positive down. 

A. ALT SEL FLIGHT TESTS 

The ALT SEL flight tests presented use the following options to invest- 
igate performance 

• Design gain set 

• Low-pitch-weight gain set 

• 30.5 m (100 ft) descent-constant throttle 

• 152.4 m (500 ft) ascent - manual throttle 

• 152-4 m (500 ft) descent - manual throttle 

• 152.4 m (500 ft) ascent - constant throttle 

• 152.4 m (500 ft) descent - constant throttle 

• Heavy turbulence 

• Moderate turbulence 

• Light turbulence 

• HDG SEL operational 

Standard operating procedure with the ALT SEL autopilot requires the pilot to 
manually adjust throttle to maintain airspeed. Flight tests with constant 
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throttle are used to test the autopilot at varying flight conditions. 

The pilot opinions for different ALT SEL gain sets are made by perform- 
ing a series of 30.5 m (100 ft) descents. Figures 13, 18, 19, 27 and 28 show 
the effect of the gain sets and turbulence conditions. Figures 27 and 28 are 
conducted with the same low pitch weight gain set but have heavy and medium 
turbulence conditions, respectively. Turbulence increases altitude hunting 
amplitudes as expected. The pilots judged the ride quality in Fig. 27 typical 
for heavy turbulence conditions. The variable, ZF, shown in the figures is 
the best indication of actual vehicle altitude estimated from a complementary 
filter using radar position and onboard accelerometers. HEIGHT ESTIMATE in 
Fig. 13 is the output of the barometric altimeter complementary filter (z(+)) 
used by the control law and HBARM is the barometric altimeter sensor output. 

Figures 20 and 31 show ± 152.4 m (500 ft) altitude changes for the design 
gain set. In Fig. 30, the autopilot is attempting to execute a 152.4 m ascent 
without increasing throttle. Airspeed steadily decreases and pitch angle and 
angle of attack increase until the pilot disengages the autopilot near stall at 
the onset of buffeting. The HDG SEL autopilot is operational during the flight. 
The autopilot did not cause the aircraft to depart from flight at a high angle 
of attack. The autopilot properly ascends when the pilot adjusts throttle as 
shown in Fig. 29. The autopilot properly descends with or without throttle ac- 
tivity. The ALT SEL autopilot has a favorable comparison to Fig. 32 which 
shows a manual 152.4 m altitude descent with constant throttle. 

The data in Figs. 28 and 29 exhibit only a slight difference in perfor- 
mance between gain sets while the pilot has a stron preference to the perfor- 
mance in Fig. 29. Evaluating performance strictly from data (or nonlinear 
simulations as is typically done) is helpful but usually does not replace 
actual pilot opinions. 

B. HDG SEL FLIGHT TESTS 

The HDG SEL flight tests presented use the following options to investi- 
gate performance: 



• 0.1 sec time constant a prefilter 



• Large airspeed variations 

• Heavy Turbulence 

• Moderate Turbulence 

• Light Turbulence 

ft Pilot entered disturbances through rudder cable 
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Multiple turns 


• ALT SEL is operational 

The investigation of different prefilters is discussed in Chapter IV, Section 
A. The flight tests results are shown in Figs. 15 to 17. The performance of 
the HDG SEL autopilot under moderate turbulence conditions with the 0.1 sec 
time constant a^. prefilter, (the nominal design) , compares favorably to the 
linear simulation in Fig. 21. The autopilot banks the aircraft, flys a co- 
ordinated turn with little sideslip in the shortest direction to the new 
heading, then rolls out of the turn just before reaching the new heading for 
a smooth capture. The HDG SEL autopilot uses the value of rudder at control 
activation as the nominal value. The value should change with airspeed since 
the nonzero rudder position is used to offset propeller effects. A flight 
test to gauge the effect of not changing the commanded rudder position is 
shown in Fig. 33. 

In Fig. 33, the HDG SEL autopilot is activated at 53 m/s in heavy turbu- 
lence. The airspeed is decreased to 35 m/s by adjusting throttle and a 45 
degree turn is executed. The HDG SEL autopilot performed the maneuver with 
no unusual difficulty due to either turbulence or rudder position and the 
pilot gave a favorable opinion of the maneuver. 

Commands to the ALT SEL and HDG SEL autopilots can be made simultaneously. 
Figure 14 shows the aircraft performing a spiral descent. Multiple spirals 
in Fig. 14 show the autopilot transferring between straight and level flight 
and turning, descending flight during the same engagement period. The pass 
through a 180 deg course in Fig. 14 demonstrates that the yaw complementary 
filter adequately suppresses the heading gyro discontinuity at 180 degrees. 
Flight tests on a 180 deg course (not presented) could not be distinguished 
from any other course heading. 

C. PITCH SEL AND ROLL SEL FLIGHT TESTS 

The PITCH SEL and ROLL SEL autopilots are simplier versions of .the ALT 
SEL and HDG SEL autopilots. The feedback gains are similiar, but not identical, 
and the command model is a simple rate limit of the amount the angle commands 
can change in 0.1 seconds. A pitch angle command of +5 degrees above the trim 
pitch angle, with constant throttle, at 18 seconds into the test is shown in 
Fig. 34. The altitude hunting in Fig. 34 compares favorably to the ALT SEL 
response during the early pitch hold phase of the flight tests. 

A roll angle command of +5 degrees is shown in Fig. 35. Prior to 
engaging the autopilot, the aircraft was trimmed with a 4.0 degree roll angle 
for straight and level flight. The ROLL SEL autopilot uses an internal roll 
command of 0.0 degrees at engage time. At ROLL SEL autopilot engage in Fig. 

35, the 0.0 roll command causes the aircraft to initiate a turn. The 5 deg 
command is initiated 28 seconds into the flight test. The roll command to 
rudder position crossfeed gain is used and the aircraft performs a coordinated 
turn. The sideslip offset shown in Fig. 35 and evident in other figures is 
believed to be sensor bias. The steady state sideslip value remains the same 

30 



before and after engaging the autopilots - 


D. APPROACH MODE FLIGHT TESTS 

The approach modes, APR GS, APR LOCI, APR LOCR, and APR LOCP are the 
most complicated of the autopilots designed in this report. An effort was 
made to take advantage of the command model formulation and design approach 
modes which work best with a Microwave Landing System (MLS) but are operat- 
ional with an Instrument Landing System (ILS) . The flight test results 
presented in this section use the Wallops Island Radar to emulate ILS signals 
only. In the future, airports will be converting to the microwave landing 
system. The MLS has a large volumetric coverage and supplies azimuth, ele- 
vation and range information to approaching aircraft. Knowing the geometry 
of the MLS transmitters, the aircraft can compute its position in cartesian 
coordinates with high accuracy, allowing the aircraft to make steep approaches 
along curved flight paths in the terminal area. In the interim period, air- 
craft autopilots will utilize either MLS or ILS. 

The ILS, which is in use today, has limited coverage with information 
supplied to the aircraft in the airspace surrounding a particular glideslope. 
Accurate deviations from the beam center line (localizer) extend 2.5 degrees 
(typically) to either side and accurate deviations from the glideslope extend 
0.7 degrees (typically) up and down. The only range information occurs when 
the aircraft receives a signal from the outer marker (about 9700 m) and from 
the middle marker (about 800 m) enroute to touchdown as the aircraft passes 
over either marker. Pilot handbooks supply information about runway heading, 
glideslope angle (2.5 - 3.0 deg.), transmitter location and actual outer 
marker and middle marker distances. 

Typical operation of a commercial autopilot requires the pilot to set 
up an intercept angle (between 45 and 90 degrees) , then arm the approach 
modes. The autopilot captures the glideslope and localizer and feeds back the 
localizer/glideslope errors. The feedback gains for the errors are some- 
times reduced as the aircraft approaches the runway to counteract increased 
angular deviation sensitivity. 

The approach modes for PIF are an attempt at designing autopilots which 
accommodate either ILS or MLS transmitted information. The approach modes 
command model use a range estimate, localizer deviation signals and glide- 
slope deviation signals to generate vertical and horizontal paths, (d and y) , 
the aircraft should follow to touchdown. The path and path errors are 
trivial to construct using MLS but are more difficult to compute using ILS. 

The primary issue is whether the range estimate shown in Eqs. 178 and 180 is 
sufficiently accurate to allow path construction. Feeding back path error 
rather than angular deviation error allows the PIF control gains to be constant 
throughout the approach mode flight. 

The operating procedure for the approach autopilots requires the pilot 
to maneuver the aircraft to the desired altitude and intercept angle for 
ILS capture. The autopilot is engaged and automatically activates the HDG SEL 
and ALT SEL modes. After the aircraft enters the localizer signal coverage, 
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the software switches automatically from HDG SEL to APR LOG. The APR LOG 
mode initializes a y-position and generates a path on line for y to approach 
0.0 using bank angle and yaw angle maneuvers and simple kinematics (Appendix 
G) . The errors between the commmand model path and y estimated from Eq. 200 
are fed back and integrated along with roll and yaw. If the intercept angle 
is large enough (greater than 30 deg) , the aircraft initially rolls away 
from the beam centerline to intercept with a tangential trajectory. 

The APR LOG mode tested is designed to capture two seconds after re- 
ceiving a valid localizer signal. If the aircraft remains banked for 
shallow intercept angles the capture sequence willl undershoot the beam 
centerline. The command model control system does not allow the undershoot 
to occur by stopping the turn sequence when a 30 degree intercept is reached 
and holding the 30 degree intercept until the capture can be continued for a 
smooth intercept. Pilot comments indicated that the autopilot should delay 
capture until this sequence can be performed with a single banked turn. This 
approach is discussed further in Ghapter VI. 

The APR GS mode command model is activated when the glideslope signal 
coverage is entered. If the aircraft is below the glideslope and the APR GS 
model control is commanding a positive flight path angle, (pitch up), the 
aircraft remains in the ALT SEL mode. When the APR GS model commands a 
negative flight path angle, the current perpendicular position between the 
aircraft and glideslope (d position) is initialized in the model and the 
model generates a path on line for d to approach 0.0 The d position is 
reduced to zero gradually so that the aircraft approaches the glideslope 
tangentially (no overshoot) . 

The flight test shown in Fig. 36 has the pilot nulling the localizer 
signal and the autopilot flying only the APR GS/ALT SEL modes. The range 
estimate is constant until the aircraft approaches the outer marker then 
propagated open loop using Eq. 178 only. This method for estimating range 
is not accurate enough after capture and the aircraft began "porpoising" 
at a range near 1600 m (1 mile) before touchdown. The pilot does not 
lower flaps during the capture and approach sequence. The aircraft smoothly 
captured the glideslope without overshoot. Altitude hunting is evident 
throughout the data in the sequence but is not noted by the pilot. 

The flight test shown in Fig. 37 uses the range estimate consisting of 
Eqs. 180 and 178. The aircraft began to smoothly capture the glideslope at 
190 seconds. As shown in the figure, the pilot lowered the flaps in two 
steps, one near 120 seconds and the other near 200 seconds. The increased 
lift due to flaps caused the aircraft to ascend. The autopilot counteracts 
the flap movement and pitched the aircraft down further at each flap change. 
The autopilot accurately followed the glideslope after the second flap 
change until the pilot took over at 60 m (200 ft.) altitude. 

The flight test in Fig. 38 is performed with the pilot flying an 
approach mode where the pilot purposely intercepted the beam centerline 
after the outer marker. The outer marker signal was not received and the 
range estimate was held at a constant value until the glideslope signal 
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was within range and both Eqs. 178 and 180 were activated. The adaptive comple- 
mentary gain in the update equation for range is used to keep large initial 
range errors from disturbing the control system. The APR GS (and APR LOCR) auto- 
pilot performed as designed even though the initial range estimate was poor. 

The pilot lowered flaps in one step coincidently with the autopilot capturing, 
the glideslope as shown in Fig. 38. 

The flight tests for the APR LOCI, APR LOCR, and APR LOCP autopilots are 
shown in Figs. 39 to 41 for a 45 degree (approx image ly) intercept. While the 
pilot favored the APR LOCI mode capture, a wings level flight along the beam 
centerline indicated APR LOCI could not accommodate the 4.0 m/s (8 kt) cross- 
wind present during flight tests. Comparing sideslip (BETA) in Figs. 39 and 
41, the APR LOCI mode sideslip gradually departs from coordinated flight as the 
aircraft tracks the beam centerline and is the source of pilot discomfort. The 
sideslip during capture indicated the APR LOCR has the worst capture performance 
while APR LOCI and APR LOCP are comparable. The APR LOCR and APR LOCP auto- 
pilots automatically accommodate the crosswind while remaining on the runway 
centerline which keeps the sideslip angle small. This is achieved without 
estimating wind magnitude or using it in the control system. The pilot could 
not distinguish between APR LOCR and APR LOCP performance and choose APR LOCR 
as the approach mode for continued flight testing. APR LOCR is similar in 
integrator command error construction to a method suggested in -Ref. 18. The 
flight test data suggest that APR LOCP is the better design. 

The APR GS flight test shown in Fig. 28 had the safety pilot track the 
localizer deviation signal. The time history of the pilot's localizer tracking 
performance is shown in Fig. 42. The autopilot's significantly greater accu- 
racy in localizer tracking with less roll angle activity is evident when Fig. 

42 is compared to Figs. 40 to 41. 

Flight test results for a 60 deg. intercept angle using the APR LOCR auto- 
pilot is shown in Fig. 43. The aircraft remains banked until intercepting the 
beam centerline then tracked with little overshoot. The pilot favored the 60 
deg. intercept capture sequence, becasue only one roll command is executed in 
this sequence. 

The localizer capture and hold sequence where the aircraft purposely 
misses the outer marker, is shown in Fig. 44. The results shown in Fig. 44 
occur concurrently with the APR GS flight test results shown in Fig. 38. The 
initial y position error is estimated to be much larger than it actually is, 
resulting in a large commanded bank angle. The bank angle oscillates as the 
changes in heading decrease the y position more quickly than expected. The 
range estimate becomes accurate 90 seconds into the flight and the aircraft 
accurately tracks the localizer error thereafter. 

The APR GS and APR LOG autopilots performed exactly as designed. Perfor- 
mance of the APR GS is excellent while the APR LOG command model requires ad- 
ditional sophistication to achieve a desirable rating from the pilot at any 
intercept angle. Favorable comparisons can be made with the autopilot ap- 
proaches and piloted approaches presented in Ref. 2 for the same NAVION air- 
craft flying the same approach procedure during display evaluation. 
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VI . CONCLUDING REMARKS 


This report provides a complete and rare practical demonstration of the 
advantages of a particular modern optimal digital flight control system design 
for use in future General Aviation autopilots. A digital control law for per- 
forming general aviation autopilot functions is designed using Linear Quadratic 
Regulator theory and Feedforward Control theory. A discrete nonlinear command 
model that accepts pilot inputs is constructed for each autopilot function and 
outputs desired aircraft state trajectories for the command states. The output 
of the command model is fed forward to a proportional-integral-filter (PIF) 
control law. The digital PIF control law forms command errors for feedback 
stabilization and output command error integration for command following. The 
PIF control law uses standard general aviation sensors for feedback and comple- 
mentary filters to smooth the output of the barometric altimeter and yaw gyro. 

A computer program is available to rapidly and completely design the PIF con- 
trol law by the general user. 

All the autopilots have been successfully programmed in FORTRAN in the on- 
board ROLM flight computer. All the autopilots have been successfully flight 
tested and performed as designed. The constant feedback gain autopilots 
stabilized the aircraft throughout the normal flight regime of the aircraft. 

The full state feedback PIF LQR design tends to be relatively insensitive to 
modeling errors. 

A. RECOMMENDATIONS 


Areas of possible improvement in the PIF design approach are discussed in 
this section and are recommended for further investigation. 

• Output Feedback - General aviation autopilots do not use full state 
feedback for stability as PIF does. General avaition autopilots use 
output feedback, which eliminates costly sensors, improves reliability 
since there are fewer sensors to fail, does not require accurate air- 
craft models to estimate unmeasured states and is easier to design 
(from a classical point of view) since there are fewer gains to manually 
adjust. Recently, in Refs. 9 and 31, numerically solutions for deter- 
mining output feedback digital control laws have been developed. The 
results in Ref. 9 are derived for the discrete LQG optimization problem 
that includes the PIF approach. The PIF control law, as designed, 
theoretically requires a full order Kalman filter for optimvim stochastic 
performance. The output feedback optimal control problem solved in Ref. 
9 formulates the LQG PIF design problem exactly as it implemented, with- 
out the Kalman filter. The resulting optimization solution would re- 
duce PIF gains for noisy sensors (ay, a^, and hg) and increase PIF gains 
optimally for accurate sensors (0, (p and r) . Complementary filters, 
actuator dynamics and sensor dynamics can be included in the PIF control 
design model. The PIF control law using an output feedback approach 
would be designed and evaluated using the implementation structure that 
exists in flight. The output feedback approach would replace ^ in Fig. 
7 with the available sensor measurements. The feedback gain C 5 and 
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the feedf 02 Tward gain in Fig. 7 would change dimension to account for 
the fewer measurements. The rest of the structure in Fig. 7 would remain 
the same. 

• Quadratic Weight Section - A new procedure for choosing the quadratic 
weights in the optimal control problem has been developed in Ref. 8. 
Designing PIF control laws using this approach should be investigated. 

• Fuel Optimal Trajectories - The command models can be made more complex 
to increase features. For example, the ALT SEL autopilot changes alti- 
tude based solely on achieving good ride quality. Properly including 
energy and fuel dynamics (similiar to approaches in Refs. 32 and 33) in 
the command model enables the command model control system to be design- 
ed to change altitude with reduced fuel cost for the aircraft. 

• APR LOC Improvements - The APR LOG command model control system should 
be improved to increase pilot acceptance. One method of improvement 
(Ref. 12) is to determine the maximum intercept angle versus speed for 
a no overshoot approach capture as a function of localizer transmitter 
distance and incorporate this into the APR LOC engage logic. 

• Flight Tests - The approach modes should be tested with a microwave 
landing system. The approach modes should be tested at a commerical 
airport using ILS. 

• Increase Robustness - New synthesis methods for increasing linear optimal 
regulator robustness have been developed at ICS and can be included in 
the PIF design procedure. 

• Stall Avoidance - Stall detection and avoidance logic should be coupled 
with the ALT SEL and PITCH SEL autopilots. Stall may occur when the 
pilot ascends without properly adjusting throttle. 

• Graceful Degradation - The ROLL SEL and PITCH SEL autopilots can serve 
as revert modes to be used automatically when sensors (heading gyro, 
barometric altimeter) fail or actuators (aileron, rudder) fail. The 
ALT SEL autopilot could switch to PITCH SEL as one method of avoiding 
stall. The ROLL SEL and PITCH SEL autopilots could be designed using 
output feedback to give deteriorated, but stable, performance using a 
minimum of sensors. 
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APPENDIX A 


AIRCRAFT EQUATIONS OF MOTION 


A. OVERVIEW 


This appendix reviews the nonlinear and linear perturbation equations of 
motion for a general aviation aircraft. These derivations make extensive use 
of vector-matrix differential equation (" state- space" ) notation, and the 
resulting equations are written in a form suitable for analysis using concepts 
of modern control theory. The notation is fairly standard but a review of 
the List of Symbols may be desirable. 

B. NONLINEAR DYNAMIC EQUATIONS 


The nonlinear rigid-body equations are reviewed in this section. The 
equations are developed using "flat-earth" assumptions, i.e., the effects of 
earth curvature and rotation are assumed negligible. This means that earth- 
fixed and inertial reference frames are equivalent. 

The origin of the inertial reference frame used here is located on the 
surface of the earth at the runway touchdown point. The x-axis points along 
the runway center line towards the approaching aircraft. The z-axis points 
down and the y-axis points to the right completing a right-hand coordinate 
system. 

Detailed vehicle state equations can best be expressed in body-fixed axes. 
These are the axes in which the pilot, the sensors, and the control surface 
locations are defined. Body axes are the only axes in which the moraent-of- 
inertia matrix is constant. Dynamic data collected from flight tests usually 
are expressed in body axes. 

The body-fixed axis system used in this report has its origin located at 
the body center of mass, and is fixed in orientation with respect to the 
vehicle. The body x-axis extends forward out the vehicle's nose, the y-axis 
extends out the right wing, and the z-axis extends out the bottom of the 
vehicle. The x-z plane is usually a plane of geometric symmetry, if the vehicle 
has one. For any nominal flight condition, another commonly used body-fixed 
axes system can be chosen so that the x-axis is aligned with the velocity 
vector, and the z-axis is in the body-axis plane. This set of body-fixed axes 
is referred to as the stability-axis system and is not employed in this report. 

A transformation from inertial to body axes is composed of a right-handed 
yaw through an angle l|j; then a right-handed pitch through an angle 9, and then 
a right-handed roll through an angle (j). The inertial-to-body axis transform- 
ation proceeds as follows: 
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Hg((j),e,i/;) = H^(6) 
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0 

-sinijj 

cosip 0 
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_sin9 

0 

COS0_ 

. 0 
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Two aerodynamic 
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in this 

report 

are the 

angle ol 


C72) 


and angle of sideslip, B. Angle of attack represents the aircraft body pitch 
angle above the velocity vector, and sideslip is the angle that the aircraft 
nose is yawed left of the velocity vector. 


The nonlinear equations of motion of the aircraft are determined by con- 
structing kinematic relationships between aircraft translational and angular 
position states and by applying Newton's Second Law to the dynamics of the 
vehicle. The kinematic relationships are 



= »B^B 

(73) 

^B 

— B 

(74) 


and the results from applying Newton's Second Law to translational and rotat- 
ional motions produces (Ref. 34) 

+ »E % - =B^B "5) 

The vector j^, is the position of the vehicle relative to the earth fixed co- 
ordinate system = [xg/ Ye/ Zg] . The inertial-body Euler angle vector is 

The Euler angle derivatives occur in three different 
references frames and are related to the body-axis vehicle rates as shown in 
Eq. 74 where 


L 


B 


1 0 -sinG 


0 cos(}) sin(f)cos0 
0 -sin(j) cos(})cos0 


(77) 


The applied specific forces consist of gravitional forces and contact forces 
as shown in Eq. 75. The specific contact force is broken into two components. 
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one of which is due to aerodynamic forces, Fg, and one of which is due to thrust, 
Tg. The gravity force is = [0 0 g] . uig in Eq. 75 is the cross product 

equivalent matrix for ^ and is given by 


r 0 -r ql 
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0 -p 


L-q p 0 J 


(78) 


The inertia matrix, Ig, in Eq. 76 contains all products and moments of inertia; 
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(79) 


while the contact moments consist of aerodynamic components. Mg, and thrust 
components, Equations 73 to 76 when combined fall into the general state 

equation form 


X = _f (2i> Jj) 

by defining the state vector 

T _ f T T T t1 
- [_-B -B -B -eJ 


(80) 


(81) 


and noting that the aerodynamic forces and moments are functions of states, con- 
trols, and to some extent, the state history. 

C. LINEAR DYNAMIC EQUATIONS 


The linearization procedure is performed by construction of a Taylor series 
expansion representing the nonlinear equations about some nominal trajectory: 

X = ^ + Higher Order Terms 


l(2o’-ao^ 


9f 

9x 


U=Uq 


A 

Ax + 

— 9u 


A^ + Higher Order Terms 


(82) 


21 = 2^0 


Here, the subscript "0" indicates the nominal value and the prefix "A" denotes 
a small perturbation. All except first-order terms are then neglected by 
arguing that the higher-order terms are small compared to linear terms. The 
results of this procedure are separated into a nonlinear equation describing 
the nominal trajectory (Eq. 83) and a linear equation defining the dynamics of 
the perturbations about the nominal trajectory (Eq. 84) : 
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(83) 


So 






Ax = AAx + BAu 


(84) 


The formal linearization of the aerodynamic forces and moments arid non- 
linear dyanmic kinematic equations is a lengthy but straight forward process 
which can be found in a number of references (Refs. 18, and 34) . The linear- 
ization of each nonlinear dynamic equation in Eq. 78 results in the following: 




Avj ' (AFj + ATg)/m + (H®^ + So*^B ' “Bo *^B 


^ <“b + *‘"b> - ®B^ 5b„S - ("bX^_ 


Aoj, 


(85) 

( 86 ) 

(87) 

( 88 ) 


The perturbation equations contain as yet only general terms for the pertur- 
bation aerodynamic effects. The only remaining undefined term is in Eq. 86, 
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(89) 


(90) 


The perturbation aerodynamic effects are assumed to be linear functions of 
the perturbation states, of the time derivatives of the perturbation states and 
of the perturbation values of the vehicle control variables. This viewpoint 
leads to the conventional name of aerodynamic stability derivatives for the 
coefficients that relate the pertrubation states, state rates and controls to 
the perturbation forces and moments. Assuming insignificant altitude and orient- 
ation effects on contact forces and moments, and assuming insignificant state 
derivative and angular rate effects on thrust forces and moments, the pertur- 
bation aerodynamic forces and moments are as follows: 
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( 91 ) 


(92) 


(93) 


(94) 


By using the definition of non-dimensional stability derivatives given 
in Ref. 24, it is possible to define the dimensional stability derivative 
matrices in Eqs. 91 to 94 in teions of the usual body axis non-dimensional 
stability derivatives. This is done in Eqs. 95 to 102. 
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The Avt^ term in Eqs. 98 and 102 is included because of the affect of d on AMpj. 

The subroutine FGAERO in the program PIFGCT, Ref. 5, accepts the non- 
dimensional stability derivatives and trim conditions shown in Eqs. 95 to 102 
and uses the information and the derivation in this section to construct A and 
B in Eq. 84. Further details of the linear model construction are given in 
documentation of FGAERO. 

Many elements in the stability matrices are known to be zero or can be 
neglected in all cases of interest. Considerable simplification of the equa- 
tions of motion for a specific vehicle can be performed and this is done for 
the NAVION aircraft in Appendix B. 
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APPENDIX B 


AIRCRAFT AER0DYN2\MIC MODEL 


The aerodynamic and mass data and simplified equations of motion used in 
the NAVION model are discussed in this appendix. Comparisons of the model 
response to the actual aircraft response for a pilot induced pulse in the air- 
craft control surface is also shovm. 

A. AERODYNAMIC AND MASS DATA 

The aerodynamic data for the NAVION aircraft is taken from Ref. 30. 

Reference 30 uses flight data rather than wind tunnel data to identify the 

aerodynamic coefficients and therefore captures rate effects due to q, r, p, 

and a,. A comparison of flight test results with linear model simulations 

caused some of the aerodynamic coefficients from Ref. 30 to be adjusted to 

better represent the aircraft dynamics. The original mismatch between the 

model and the aircraft was sufficient to cause the lateral-axis PIF control 

law to be unstable. The aircraft wheels in current flight tests are in the 

takeoff and landing configuration (increased drag) while the wheels in Ref. 30 

are believed to be retracted. The NAVION airplane parameters are shown in 

Table 14. The gross mass value represents the best estimate of the current 

fully instrumented NAVION. The mass value (WT) is 204.8 kg (452 lbs) heavier 

due to instruments than the value given in Ref. 30. The instruments are near 

the c.g. The values for I , I , and I remain unchanged from those given in 
- ^ ^ X Y Z 


The matrix differential equations for the nonlinear aircraft model in 
Appendix A can be simplified to the following series of scalar nonlinear 
equations for the NAVION aircraft. 


u = rv - qw - sin(0)g + % pV^A (103) 

V = pw - ru - cos(0) sin((f))g + pV^A (104) 

w = qu - pv + cos(0) cos((}))g + ^ pV A (105) 

p = “ ((I -I^)qr + pV^A b C^) (106) 

q = ^ ((I -I )rp + % pV^A c C^) (107) 

^ I z X ^ ^ M 

y 

f = ^ (d^-I )pq + is pV^A b C^) (108) 

z ^ 
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0 = cos(<j))q - sin((J))r 


0 = p + (q sin((J)) 4- r cos((j))) tan(0) 


ip = (sin((p)q + cos(<J>)r)/cos(6) 
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(109) 

( 110 ) 
( 111 ) 
( 112 ) 

(113) 

(114) 

(115) 

(116) 
(117) 


Two flight conditions are used in Ref. 30 to extract aerodynamic parameters. 
They are shown in Table 15. The wheels down fully weighted NAVION used in 
flight tests cannot fly as fast as condition I- Corrections based on flight 
tests data have only been performed for condition II. The longitudinal aero- 
dynamic parameters at the two flight conditions are shown in Table 16. The 
lateral aerodynamic parameters at the two flight conditions are shown in Table 
17, The two parameters which have been adjusted after flight tests are Cy^ 
(from -0.68 to -0.143) and C^g^ (from -0.007 to -0.023). ^ 

Reference 35 was obtained after the completion of the PIF design and 
flight testing. The analog match from Ref. 35;, and shown in Tables 16 and 17, 


44 



is for the second NAVI ON, a variable stability aircraft equipped with side force 
pannels on the wings. Cyg and from Ref. 35 are in general agreement with 

the modifications used in the PIF design. The analog match in Ref. 35 did not 
use Czq and as identifiable parameters and for this reason, the other 
longitudinal parameters in Table 16 do not agree well with Ref. 30. 

B. SIGN CONVENTION 

In implementing a control law it is crucial that the sign conventions 
used during the control law design be explicitly detailed. For the following 
presentation of sign conventions, the aircraft is assimed to be near the 
glidepath approaching the runway. Left and right is referenced to the pilot 
perspective while setting in the cockpit. Vertical height and normal acceler- 
ation plotting from flight test are opposite in sign to the convention used in 
this section. The barometric altimeter complimentary filter also has the 
opposite sign convention. 

The earth fixed (inertial) coordinate system and body axis coordinate 
system sign conventions are discussed in Appendix A. The following are sign 
conventions for positive step changes in control variables. The sign con- 
vention is chosen to match the sign convention onboard the NAVION aircraft 

Positive Rudder - Positive Yaw 

Positive Yaw - Aircraft Nose Rotates Right 

Positive Rudder - Rudder Trailing Edge Rotates Towards Right 

Wing 

Positive Aileron - Positive Roll 

Positive Roll - Right Wing Down 

Positive Aileron - Left Aileron Trailing Edge Rotates Down, Right 

Aileron Trailing Edge Rotates Up 

Positive Elevator - Positive Pitch 

Positive Pitch - Nose Up 

Positive Elevator - Elevator Trailing Edge Rotates Up 

As the aircraft progresses down the glideslope the flight path angle, y 
is given by 

Y = sin"^ (-z/V) ( 118 ) 
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V is the aircraft total earth relative velocity and z is the aircraft's earth 
relative vertical velocity. If the aircraft is descending the flight path angle 
is negative- When the aircraft has no lateral error, the distance in the verti- 
cal plane from the aircraft c.g. in a perpendicular line to the glidepath is d, 
and satisfies the equation 


d = V sin(X^g + X) 


[119] 


Figure 45 shows the aircraft below the glideslope. The variable Xqq is the 
glideslope angle and is defined to be positive for typical glideslopes. Hence 
for a desired glideslope angle of 3 deg, the aircraft's flight path angle must 
be -3 deg for the aircraft to be tracking the glideslope, i.e., d = 0. If the 
aircraft is above the glidepath, d is positive. The angular error, F, relative 


to the glidepath is 
r = sin ^(d/R) 


where R is the range to the aircraft, 
r is positive. In an MLS environment, 
using the equation 

d = X sin(Xgg) - z cos(A^g) 


(120) 


If the aircraft is above the glidepath, 
d can be computed from x and z position 


( 121 ) 


The localizer beam centerline intercept and hold geometry is shown in Fig. 
46. Projecting the velocity vector onto y produces 


y = V cosy sin(?-^^^^) 


(122) 


where ^ is the velocity heading angle measured relative to a reference value 
?ref* lateral distance y satisfies the equation 


y = sin(H) 


(123) 


The localizer deviation angle, 5, and the distance y shown in Fig. 46 are 
negative. The distance is the projection of R onto the horizontal plane. 

If the assumptions are made that the flight path angle and sideslip angle 
are small, (and the aircraft is in the approach mode) then the above equations 
can be simplified to 
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y - V sin(Tp-i/j^^^) 


(124) 



I 


y - R sin(H) [125] 

The lateral position, y, and rate, y, can be estimated using Eqs. 124 and 125 
and sensor information available onboard the aircraft. Projecting the velocity 
vector onto R and using the same simplifying ass\jmptions produces 

R - V - E) [126] 

C. SIMULATION AND FLIGHT TEST COMPARISON 

During early flight tests, the safety pilot purposely made separate small 
pulse inputs to each aircraft surface. The resulting NAVION aircraft response 
from the sensors was recorded for comparison with the linear open- loop model 
simulation response for similiar pulse inputs. The results of the comparison 
for an elevator input and a rudder input are shown in Figs. 4 and 3, respec- 
tively. 

The normal accelerometer sampled output (10 samples per sec) and pitch 
rate gyro sampled output for an elevator input agree well with the linear simu- 
lation as shown in Fig. 4. The normal accelerometer output sign convention is 
positive up and biased by -1.0 g to give zero output in level flight. 

The aircraft response to a rudder pulse originally had poor agreement with 
the model response. After adjusting two aerodynamic coefficients as discussed 
in Section B, the new comparison in Fig. 3 is considered to be adequate. The 
initial response for the first few seconds match well then begin to diverge as 
the simulation continues. The cross-axis aircraft response of roll rate for a 
rudder change has less agreement with the simulation than sideslip and the 
lateral accelerometer output. The aircraft's Dutch Roll mode appears to be 
better damped than the model response. The small rudder movements made by the 
pilot during the pulse command may have contributed towards the simulation and 
aircraft disagreements. The linear model simulation disagreements did not ad- 
versely affect the PIF control law performance in flight. On the other hand, 
better identification of all aerodynamic coefficients should further improve 
PIF control law perfo 2 miance. 

D. SENSOR MEASUREMENT PERTURBATION EXPRESSIONS 


The PIF control law employs, as much as possible, sampled sensor outputs 
for direct feedback. The output of some sensors, particularly the accelero- 
meters, the pitot tube, the sideslip vane, and the angle-of-attack vane, are 
nonlinear functions of the aircraft body-axis states used in PIF feedback. The 
PIF design procedure at a trim flight condition requires an expression for the 
perturbation sensor output in the form 


Az = H Ax + D Au 


(127) 


where H and D are constant aircraft state control observation matrices. The 
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purpose of this section is to derive perturbation expressions in the form of 
Eq. 127 for the previously mentioned sensors. Sensor inaccuracies and noise 
effects are neglected. 

The output of an orthogonal traid of body-mounted accelerometers aligned 
with body coordinates is expressed as (Ref. 3) . 

Sb ' is “b + [“b “b ^bJ % (128) 

A flat, non-rotating earth surface is assumed. The gravity vector, gjg, the 
Euler angle transformation matrix H®, the cross product notation UJg, and the 
body-axis velocity vector, are defined in Appendix A. Axr is the position 
of the accelerometers with respect to the center of gravity in body axis. 




Ax 


Ay. 


Az 


(129) 


Perturbation expressions from individual componets in Eq. 127 are derived in 
Ref. 3 except for the third term which follows. 


A I jj® u>l + ajjAxg) - a(-m® Aig (jj®) + a(-A:j li® 
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The perturbation expression for the accelerometers becomes 
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(131) 


where Lg is defined in Eq. 77. The variables Avg and A^ are part of the 
perturbation aircraft dynamics. 


Ax = A Ax + B Au 


(132) 


discussed in Appendix A. 

Partitioning Ax, A and B, it follows that 
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(133) 


Av^ = A Av + A Ao)_ + A ,Av^ + B Au 
— B w ~B VO) — B vv — B vu — 


Aw^ = A Av^ + A Ati)^ + A Av_ + B Au 
— B lov — B 0X0 — B 0 )V — B 0)u — 


(134) 


Substituting Eqs. 133 and 134 into Eq. 131 and combining terms results in the 
perturbation measurement equation for the body mounted accelerometers 
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[135] 


If Eq. 135 is used as the measurement equation in constructing the PIE gain 
matrices then a sampled accelerometer output can be used directly as a feedback 
state in the control law without adding in gravity or correcting for off center 
of gravity locations. There are certain advantages which can be obtained using 
off c.g. located accelerometers as discussed in Ref. 29. Although it was not 
feasible, the problem discussed in Chapter IV, Section C could have been 

alleviated by moving the lateral accelerometer forward along the body x-axis. 
Corrections for gravity and off c.g. locations must be performed when the ac- 
celerometer outputs are used in complementary filtering as discussed in Chapter 
IV, Section A. 


The pitot tube measurement of velocity magnitude, V, the sideslip vane 
measurement, B, and the angle of attach vane measurement, a, are related to the 
body-axis velocities by 


V 


■^u2 + v2 + 

B 

= 

tan~-^(v/ju2 + w^ ) 

a 


tan“l (w/u) 


Perturbations of the above equation results in 


AV 

AB 

_ Aa. 






"Au“ 

Av 

_ Aw_ 


[137] 


Jw is a diagonal matrix which has elements 1.0, V^, and V^ cos Bq* Hg is the 
body to wind axis transformation and is the transpose of 
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cosa 0 -sina cos3 sinS 0 
o o o o 


0 1 


sina 0 cosa 
o o 


-sin3 cos6 0 
o o 


Equation 137 can be used to form H and D in Eq. 127 when the air data sensor 
information is needed for feedback. 
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APPENDIX C 


AUTOPILOT MODE MODELS AND CONTROL SYSTEMS 


Each autopilot mode is composed of a dynamic model and a dynamic model 
control system which are described in this appendix. When an autopilot mode is 
engaged, the internal states in the command model are initialized to the air- 
craft states. After initialization, the command model is propagated to generate 
state trajectories the PIF control law follows. Some variables in the command 
models (airspeed, range, pitch angle) are continually updated to match aircraft 
conditions. Other variables (yaw, roll, height, lateral position) are internal 
to the command model. 


When the ALT SEL, HDG SEL, APR GS and APR LOG command models receive a 
pilot command, a nonlinear feedback control system maneuvers the command model 
to the new conditions. The feedback control system allows the autopilot command 
models to operate over a wide range of initial (boundary) conditions in real 
time. The nonlinear feedback control systems are designed solely to provide 
good ride quality. The simple autopilots BETA HOLD, ROLL SEL and PITCH SEL 
operate open loop. The PIF design process requires that linear representations 
of the command model be used, hence, both linear and nonlinear models, where 
applicable, are given in each of the following sections. The basic linear com- 
mand model representation is shown in Eqs. 6 and 7. 

A. 3 m, 4)m COMMAND SYSTEM - BETA HOLD 


The sideslip, roll angle, command system is used to compute the 

roll command to rudder command crossfeed gain used in ROLL SEL. The crossfeed 
gain, aj^, is computed as an element in the feedforward matrix A 22 9. 

The linear and nonlinear command models are identical for this command system. 


^,k+l 


0 X . + 0 

— m, k 


m 
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^,k 


0 X , + I 

— m,k 


m 


m 


BETA HOLD is not intended to be an inflight command system. 


(139) 


[140] 


^m^ *^rm COMMAND SYSTEM - ROLL SEL 


The roll angle, ())^, rudder position, 5^^^, command system performs the 
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standard wing leveler plus yaw damper autopilot feature. The linear and non- 
linear command models are identical for this command system. 


•^,k+l 


0 X , + I 
— m,k 



(141) 



The crossfeed gain, aj^, is scheduled as a function of airspeed as shown in Eq. 
70. The PIF autopilot with the ROLL SEL command model is designed to accept a 
roll command, (usually 0.0) and a rudder command, <Sj^, (usually rudder pedal 

trim position) and hold these values. The autopilot will also capture and hold 
a nonzero roll command input from the pilot. 

C. 4> + Ip , 5 COMMAND SYSTEM - HDG SEL 

m ^m rm 


The Heading Select and Hold command system performs the bank-to-turn auto- 
pilot feature. When the autopilot is turned on, it holds the current aircraft 
heading of the aircraft until a different heading is requested. When a new 
heading is requested, the autopilot banks the aircraft into a steady coordinated 
turn, then establishes and holds the new heading. 

The bank to steady turn feature of the autopilot is accomplished using the 
basic nonlinear kinematics 


[b = lb 

m,k+l m,k 


At tan (() , 


[143] 


d) = d) 

m,k+l m,k 


+ At (j) 


m,k 


[144] 


The nonlinear relationship between derived in Ref. 34 and attempts 

to produce a coordinated turn where the areodynamic force lies in the body-axis 
xz plane. The difference between actual and command heading is used to form a 
roll angle command 


Zi, - 
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The command error resembles a typical feature of a nested autopilot design where 
the heading error is used to compute a roll command in the ROLL SEL mode to pro- 
vide a heading select capability (Ref. 18) . 


The linear perturbation command used in the PIF design model is obtained 
by linearizing Eq. 143 and regrouping Eq. 145. The roll dynamics in Eq. 144, 
as discussed in Chapter 3, Section A, are not represented in the linear model 
because poor transient behavior occurs if is used as the linear models con- 
trol input, Uj^. The linear model is 


Aijj 
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(148) 


The HDG SEL nonlinear command model uses ^ model control variable 

to cause 'IJitirk to smoothly transfer from one pilot requesting heading command, 
to another. The nonlinear command model control system used to follow the 
pilot command inputs is as follows: 


IF FIRST PASS THEN , 

m,o k 

m,o k rm,o 

= 'r,k 


’ '<'c,k - ’f'm.k 

; compute command 

model yaw error 

(149) 

uc = 

; compute command 
control command 

model feedback 

(150) 

4> = uc 

m,k 

; initialize roll 

rate command 

(151) 

IF Id) l>d) THEN i , =0.0 

' m,k‘ ^max m,k 

; set roll rate to zero if in 
steady turn 



; reset roll rate command if 
approaching desired heading 


^max2 = 



(152) 
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IF Ip^ AND ^ , ARE OPPOSITE SIGN THEN i o J v 

S m f ic ^ xti f ic 

• • • 

(b = MAXIMUM(cj) <J> ,) ;apply full rate command if (153) 

^max max2 maxi 

opposite sign occurs otherwise 
hold roll rate command below 
^maxl 

IF then ,,=SlGN(i 

in f iria.^ 

; check to see if low velocity (154) 
is decreasing turn radius 

;keep roll angle below ({>^^0x2 (155) 

unless <j)maxl indicates it 
should be lower 


The control system feeds back the model yaw error and model roll angle using the 
gains HSKj^ and HSK 2 . If the model roll angle exceeds a specified value/ the con- 
trol/ (}’in,k^ I® reset to zero. During this phase the model dynamics are in a 
steady turn. If the roll rate control command causes the yaw angle to accelerate 
in the direction of the yaw angle command/ ip^^, the roll rate control command is 
limited to the value 4>inaxl* angle error is below a threshold/ ^cloSE' 

full yaw angle deceleration through the roll rate control command is used. The 
last feature causes the model and the aircraft to anticipate the roll back to 
zero just before completing the turn so that the yaw angle overshoot is small. 

A minimum radius of turn/ Riji/ is used to reduce the maximum allowed back angle/ 
^maxl' filtered airspeed/ V, drops to low values. The computed value of 

^m,k' after Eq. 153/ is used in Eq. 144 to generate the trajectory. The 
resulting 4>ni trajectory is used in Eq. 143 to generate the trajectory in 
flight. 

D. z COMMAND SYSTEM - ALT SEL 
m 


m , X max 


m/X 


(p T = tan ^(V.^/R_*g) 
^raaxl ' k T ^ 


Cj) = MINIMUM ((}) , , ^) 

max maxi max2 


The Altitude Select and Hold command system allows the pilot to hold the 
current altitude or select a new altitude the aircraft should establish. When 
the pilot selects a new altitude/ the autopilot slowly accelerates (or deceler- 
ates) the aircraft into a steady climb (or steady descent) until the new alti- 
tude is reached then captures and holds the aircraft at the new selected altitude. 

The altitude select feature of the autopilot is accomplished using the 
basic command model 




"1.0 

At ' 


1 

N 

3 



+ 

> 

0 

0 

1 



_0.0 

0.1, 


^m,k 


L At J 


— k -^/k ^k ^m,k 


[156] 

[157] 
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The altitude acceleration, Zni,k' used as a control variable to transfer 
from one desired altitude value to another. The linear command model used in 
the PIF design does not use Zj^, and has the following simple form 


Az , , T = Az , + At z , 
m,k+l m,k m,k 

(158) 

Aik “ 

(159) 

Aim,fc ^^m,k 

(160) 


The altitude select command model control system is as follows where z^ 
is the requested altitude made by the pilot. 


IF FIRST PASS 


THEN z = z, , z =0.0 
m,o k m,o 


= z 1 ~ z . 

E c,k m,k 


; compute command model height (161) 


UC = ASK*z-ASK*z , 
1 E 2 m,k 


; compute command model feedback (162) 
control command 


z , = UC 
m,k 


; initialize acceleration command (163) 


IF z , >z THEN z =0.0 

' m , k ' max m , k 


; set acceleration to zero if in 
steady descent (or ascent) 


IF z ^z^^ THEN z , =UC 

' e' CLOSE m,k 


; reset acceleration command if 
approaching desired height 


z ^ = 0.0 

max 2 


(164) 


IF z„ AND z , ARE OPPOSITE SIGN THEN z ^ , 1 

E m,k max2 ‘ m,k' 


Z = MAXIMUM (z z ) 

max max2 maxi 


; apply full acceleration if (165) 

opposite sign occurs otherwise 
hold acceleration command below 

^maxl 


The control system feeds back the model height error and model vertical velocity 
using the gains ASKj^ and ASK 2 . If the model vertical velocity exceeds a spe- 
cific value, the control, Zni,k' reset to zero. During this phase the model 
dyanmics are in a steady descent (or ascent) . If the vertical acceleration 
command causes the model height to accelerate in the direction of the vertical 
height command, the vertical acceleration control command is limited to the 
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value Zmaxl* height error is below a threshold, ZcLOSE' full deceler- 

ation through the vertical acceleration control command is used. The computed 
value of Ziti^k after Eq. 165 is used in Eq. 156 to generate the Zm,k trajectory 
in flight. The altitude select and heading select model control systems are 
coincidently structurally identical. 

E. 6 COMMAND SYSTEM - PITCH SEL 
m 


The pitch angle command systems holds the current pitch angle or commands 
a new one using the simple model. 


^,k+l 


0. 0 X , + 0.00 1 

— m,k m,k 


(166) 


y,=0.0x ,+1.06 1 

-^,k — m,k m,k 

(167) 

F. GLIDESLOPE CAPTURE AND HOLD - APR GS 


The Approach Glideslope autopilot mode captures and tracks the glideslope 
The switch from ALT SEL to the APR GS mode occurs automatically when the air- 
craft is approaching the runway. The ALT SEL mode must be engaged and the air 
craft must be holding a desired altitude prior to capturing the glideslope. 
Glideslope capture can be performed either above or below the glideslope. The 
localizer capture and track should occur prior to glideslope capture. These 
requirements are standard in most autopilots. 

Glideslope capture and track is accomplished 

using the nonlinear model 

' '^m.k + \ ^m,k> 

(168) 

A - ^m.k ' \ '’k - '^m.k = \ - "“m.k 

(169) 


The distance, d^^k perpendicular distance from the aircraft to the glide- 

slope, Agg is the glideslope angle, A^^^ is the commanded aircraft flight path 
angle and V is the ground speed. In the PIE control law, Eq. 168, is used to 
define a path that the measured value d tracks to zero. The values used for V, 
Aq, and Ajjj do not have to be accurate in flight since inaccuracies only imply 
a perturbation of the trajectory d^ traces to zero. In flight, V is obtained 
using measured airspeed, Agg is initialized by the pilot and Ajj, is a control 
feedforward variable. A measured aircraft flight path angle is not required 
for feedback. The linear perturbation command model used in the PIF design is 
obtained by linearizing Eq. 168, 


Ad T = Ad . + At V, AA . 

m , k+1 m , k k m , k 
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(170) 


■HU I I I 


I 






" *^m,k 


(171) 


= Adj^ (172) 

The APR GS nonlinear command model uses a control variable to 

cause djjj^k smoothly transfer to zero from the measured value of d that occurs 
when the glideslope engage logic activates the APR GS command model controller. 
The command model control system used to construct the path, which includes the 
range estimation equations, is as follows 


DATA GSENGAE/FALSE/ 
DATA UC/0.1/ 

DATA OMKRCK/FALSE/ 


; initialize APR GS flag 
; initialize US to 0.1 
; initial outer marker flag 


IF OMKRSOLID 


THEN OMKRCK=TRUE 


/change flag when aircraft passes 
over outer marker 


IF NOT OMKRCK AND NOT GSENGAE 


THEN R|^=0UTMARK /initialize distance estimate 
to outer marker distance 


ELSE 

BEGIN 

RK = 0.0 


IF LOCENGE THEN 


RMEAS = \(t)/Xgg 


IF GSENGAE THEN 


RK = r * z, + r^ 
1 k 2 


IF RK < 0.0 THEN RK=0.0 

IF RK > RKMAX THEN RK=RKMAX 
Rj^ = Rj^ + RK * (RMEAS - Rj^) 


/initialize range estimate 
filter gain to zero 


/propagate distance extimate 
if H is a valid signal 

/ estimate R using height 
estimate and glideslope angle 


/ increase RK as aircraft 
descends on glideslope 


(173) 


(174) 

(175) 

(176) 

(177) 


/update distance estimate 


(178) 


END 
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IF R, < R . 

K min 

THEN ^ =R . 

k min 

;keep distance estimate above 

R . 

min 


IF |r, 1 < r 

' k ' max 

THEN 



BEGIN 






; calculate d estimate when 
glideslope is in range 

(179) 

IF UC > 0.0 

THEN d , =d, 
m,k k 

; reset d^^^ to d until the APR GS 
model command control system 
causes a pitch down command 


IF UC < 0.0 

THEN GSENGAE= 

=TRUE 


d^ = (0.0 - d 
, E m. 

* 3.281 

; convert to feet 

(180) 

GSK = g. |d , 1 
1 ' m,k' 


; compute APR GS command model 
control system gain 

(181) 

IF GSK > g^ 

THEN GS K==g^ 



IF GSK < g^ 

THEN GSK=g, 
4 



UC = GSK * * 

0.01745 

; convert to radians 

(182) 

UC = + UC 

GS 



(183) 

IF UC ^ A 

max 

THEN UC=A 

max 

;A should be 0.0 

max 


IF UC < A . 

min 

THEN UC=A . 

min 

/•keep flight path angle command 

above A . 

min 


A = (UC - A , 
m m,k- 



(184) 

IF |A I > A 

THEN 



A = SIGN (A ) * 
m m 

A 

max 

/•keep flight path angle rate of 
change below A 

max 

(185) 

A , = A V 1 

m,k m,k-l 

A 

m 

/•propagate command model 

(186) 

d , = d , + At V, sin (A „ + 
m,k m,k k GS 

m,k 

(187) 


END 


At the beginning of the control law computation/ the glideslope angle is initial- 
ized by the pilot. The distance, R, is initialized to a typical outer marker 
distance. If the aircraft has intercepted or passed the outer marker and the 
localizer deviation angle is in range, the distance R is propagated. When the 
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aircraft is descending on the glideslope an alternative estimate of R is made 
using the barometric altimeter z estimate and the geometry involved. The pro- 
pagated R and estimated R are combined in a complementary filter fashion. The 
filter gain, RK, is increased as the aircraft descends. When the estimated R 
falls below a prespecified value, R is fixed to the value. Equation 180 forms 
the command model control error. When -GSK * d^^k less than the flight 

path angle command, X, begins to command the aircraft to pitch down even when 
the aircraft is below the glideslope. The aircraft follows a d trajectory which 
intercepts the glideslope with no overshoot. The gain GSK is programmed to 
increase as the aircraft closes in to the glideslope. The computations after 
Eq. 183 check to determine if the commanded flight path angle is within spec- 
ified limits. Setting Xj^ax zero ensures the command only executes pitch 
down glideslope captures. The rate of change of the flight path angle is also 
checked to keep the vertical accelerations small during capture. 

G. LOCALIZER CAPTURE AND TRACK - APR LOCI, APR LOCR, APR LOCP 


The localizer capture and track autopilot mode is the most complicated 
command model design in the PIF control entourage. A number of options are 
available to specify a command model which captures and tracks the localizer 
beam centerline. Each option has advantages and disadvantages which are de- 
scribed in the flight test discussion. Chapter V. Three APR LOC command modes 
are presented (APR LOCI, APR LOCR, APR LOCP). PIF control law designs for all 
three modes are discussed in Chapters IV and V. 

The APR LOC autopilot captures and tracks the inbound localizer beam 
centerline. All three approach localizer autopilot designs use the same basic 
command model. The autopilots differ in the construction of the command error 
signal which the PIF control law must regulate to zero. 

The APR LOC mode requires that the aircraft be in HDG SEL mode prior to 
engage. The engage logic for transferring from HDG SEL to APR LOC is performed 
when the localizer signal is in range. The internal intercept angle in the 
command model control system during capture is fixed. The APR LOC command 
system has considerable room for improvement as discussed in the recommendations 
section of Chapter 6. 

Localizer capture and track is accomplished using the nonlinear dynamics 
shown in Eqs. 124 and 143 where, for simplicity, the reference value is assigned 
to be 0.0, 


y 1 J.1 = y 1 + At V. * sin ip , 
•^m,k+l -^mjk k m,k 


ip T, . 1 = ^ 1 + At tan (f) , 

’^mjk+l ^m,k V-^ ^ra,k 


[188] 

[189] 


The three APR LOC modes have the following command errors 
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APR LOCI 


- ^,k == 
APR LOCR 



1.0 o.o-j 

a, 1.0 6 

k mr , k 


(190) 


(191) 


APR LOCP 

^k -^,k = 




(192) 


The roll angle command, is the control variable for the nonlinear command 

model control system. 

The APR LOCI mode forms the command errors shown in Eq. 190 using roll 
angle and y position. The aircraft's y position is estimated from H measurements 
using Eq. 123. 

The APR LOCR mode shown in Eq. 191 is simply the HDG SEL mode with an addi- 
tional modification. The ^ni,k command in Eq. 145 is replaced with ^"ky {yk-^m,k) • 
An error in y position causes a angle error which in turn causes a coordinated 
roll command to bank the aircraft in the direction to null the y position error. 

The APR LOCR mode is similiar in principle to a localizer autopilot mode dis- 
cussed in Ref. 18. 

The APR LOCP mode is similiar to the APR LOCR mode except that the error 
is removed and only y position error commands the bank angle. The state re- 
mains in the model dynamics and is used in the PIF control law to form the x-x^* 
error (see Eq. 57) for ip. The APR LOCP mode is able to perform a coordinated 
turn, intercept the localizer beam centerline and track the beam centerline 
crabbed, if necessary. 

The linear perturbation command models used in the PIF designs are obtained 
by linearizing Eqs. 188 to 192, 
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1.0 At V cos i|) * (1.0/57, 

o o 



0.0 rA 4 ) 


At(g + 4,^ tan ^5 


APR LOCI 


0.0 0.0“1 fAy 


Ll.O 0. 


oj AI|J^ 


'l.O O.olfAcj) , 

m y iC 

_ 0.0 O.oJ A6 . 

rm,l 


APR LOCK 




-k.k 

ijj y 

-a, k.k 
k y 


-k, Ay 
0.0 Aifj 


rO.O 0.0“1 [Acf) . 

m,k 

0.0 l.oj A6 

^ [_ rm. 


*'f’k ■ % *'*’k '"'I- N 

*«r.k - *k '^4. '^♦k ■" “’k % '^y *^k 


APR LOCP 




-a, k 
k y 


JUi i.. .JL: 


■ N ^^k 

*^.k - “k N ^’'k 
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The APR LOG autopilots use as a control variable to cause to smoothly 
transfer to zero from the measured value of y that occurs when the localizer 
engage logic activates the APR LOG command model controller. The command model 
control system used to construct the path is as follows: 


DATA LOGENGE/FALSE/ 


; initiaze APR LOG flag 


= (Rj^ + RLOG) sin 


; estimate the y position using (200) 
the range estimate 


IF <c 

' k ' max 


THEN LOGENGE=TRUE 


IF NOT LOGENGE THEN y , =y, 

m,k k 

d) =(b 

^m,k 

/\ 

ih =ip 

^m,k ^k 

rm,k rk 


IF LOGENGE 
BEGIN 


THEN 


; reset the internal command model states 
until 5 is in range 




/•compute APR LOG command model (201) 
control system gain 


IF LOG ^>^2 
IF L0G^<£^ 


THEN L0G^=£^ 
THEN LOC^=l^ 


LK1=L0G, (V *3.28) 
1 k 


;try to keep closed-loop system (202) 
invariant by varying control 
gains with velocity 


LK2=LOG2/(g/Vj^) 


(203) 


Ip =LK *(0.0-y ,) 

^c 1 m,k 


; compute command 


(204) 


IF 111; |>l|j THEN 

' c ' cmax 




UG=LK *0p -1|; ) 

2 o in / iC 


;ip^ should be less than the 
4’cmax iJ^tercept angle 

; compute <J>nj command (205) 
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IP luclxl) 

' ' max 

THEN UC=sign (UC) *d) 

^max 


(j) =uc-4> , , 

m m,k-l 

;keep roll and roll rate below 
maximum values 

(206) 

IF |(J) |>d) 

' m ' ^max 

THEN d) =sign (d) ) *d) 

^m ^ ^m ^max 


d) =d) +(b 

m,k m,k-l m 

;the computed roll command used 
in Eqs. 188 and 189 

(207) 


END 


The localizer model control system will begin to intercept the localizer 
beam centerline if the pilot has activated the approach mode (APR = TRUE) and a 
valid localizer signal is available. The pilot must set up an intercept angle. 
The APR LOG mode performs best if the pilot orients the aircraft at a ~60 deg 
intercept just before the outer marker. The command model values for 
and are continually updated to the corresponding measured aircraft's values 
as long as the APR LOG control system has not been activated. Equation 204 forms 
the command models heading command whose upper value is limited to the intercept 
angle set in the control system. As the aircraft nears the beam centerline the 
heading command gradually changes to the reference 0.0 deg runway heading. The 
measured heading is referenced to the runway heading. After the heading command 
is computed, the rest of the control system resembles the HDG SEL command model 
system discussed in Section C. 
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APPENDIX D 


TRIM, STEADY STATE AND THE STAR TRAJECTORY 


This appendix discusses trim for nonlinear systems, steady state for 
linear systems and the relationship between the two. Steady state for linear 
systems is generalized to include the theory of feedforward control and model 
following. Feedforward control theory is used in the PIF derivation in 
Chapter IV. 

For a constant command, 

V = h(x , u ) (208) 

o — O — o 

= 0.0 (209) 


an aircraft is in static trim if there is a u such that Eqs. 208 and 209 are 

— o 

true and the nonlinear aircraft dynamics, shown in Eq. 80, satisfy 



u ) 

— o 


(210) 


The outer-loop states x, y, z, ifj are not included in Eq. 210. For the pertur- 
bation linear system, static trim or steady state is similiarily defined 


0 = 

A Ax’*! + B Au* 

(211) 


= H Ax’*’ + D Au’*’ = A_y^ 

(212) 

A^* 

= 0.0 

(213) 


Grouping the linear equations, static trim for linear systems can be 
easily computed by solving the linear matrix equation 
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(214) 



for Ax* and Aii*. Ax* and Au* are defined to be the "star trajectories" of 
the plant for constant A^. A solution to Eq. 214 exists if the plant has 
no transmission zeroes at zero, i.e., the quad partition matrix in Eq. 214 
has full rank. Using the definition of perturbation systems, and the Taylor 
series expansion, different nonlinear trim points caused by perturbation 
changes in represented as A^, are related through the following equations. 


y + Ay = y + Ay* = y (215) 

-*^o ^ o,new 


X 

+ 

Ax* 

= X 

(216) 

^o 



— 0 ,new 


u 

+ 

Au* 

= u 

(217) 

— o 



— o , new 



Steady-state for discrete systems is similiarly derived. The state trans- 
ition equation for a linear system is. Ref. 36, 

Ax(t) = e^^^"^°^Ax(to) + / BAu(t) dT (2is) 


Starting from time t , Eq. 218 describes how Ax(t) evolves if Aji(T) is known. 
If Aja(x) is held constant over the constant interval At (a zero-order hold 
implementation) then Eq. 218 can be represented at each time instant t, , where 


^2k+l = 


(219) 


which is an exact representation. 

The discrete PIF control law uses Eq. 219 as the discrete plant model. 
In steady state, Eq. 219 can be combined with Eq. 212 to produce a result 
sirailiar to Eq. 214, 


'(<t>-i) r' 

'l' 

1 

11 

1 

o 

1 

. H dJ j 

Lau*J LiJ 


% 


( 220 ) 
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If the discrete plant has no transmission zeroes at 1.0 then Ax* and Au* 
can be computed by 


Ax^ 

- 0 ■ 

Au^ ^21 ^22 



where 



1 

^11 

^12 


■ I 0" 

H 

dJ 

-^21 

^2 2 


_0 I_ 


C221) 


( 222 ). 


The steady-state star trajectory computed for the discrete-time plant model 
is identical to the values computed using Eq. 214. 

Tracking a constant command can be generalized to command model output 
tracking. The discrete-time plant representation for the command model shown 
in Eq. 7 and 8 is determined using Eq. 218, where A;^ is constant between 
sampling intervals, and has the following form ' 


*Sm,k+l ‘ ‘I’m ^,k + 


(223) 


The solution for the control, which causes A^ to track A^^j^ for the 

discrete-time command model is accomplished by generalizing Eq. 220. The 
generalization primarily occurs by interpreting A^ and A^ as ideal traj- 
ectories (the star trajectories) that occur when Ay* = A^^. 

To derive the equation for the star trajectories, the star trajectories 
are assumed to be related to the discrete model through constant matrices 



ll’ 

21 


12 


22 


Ax 
— m 


' m, k 


1 


(224) 


The matrices A. . in Eq. 224 are the feedforward matrices and remain to be 
determined. 


The derivation for a solvable expression for the feedforward matrices 


proceeds by forming two different expressions for ^ 


■Axtn - *251 


then 



equatiflg them. One expression is obtained using the plant dynamics 


- 


Ax*. 1 
— k 



r“ 

■^11 

^12 


1 

_ 


H 

D_ 

A 

L-^21 

^22- 


1 Au , J 

^ — m > k 


( 225 ) 


The other expression is obtained using Eg. 224, 


- -! 


A..(<p -I) A-.r A 
11 11 m 


H 


m 


D 


m 


1 

12 

0 


Ax 

"m,k 

Au , 
—m,k 

Au , ^ 
— m,k+l 


Au , 
-m,k 


(226) 


To proceed with the derivation, Aj^^k +1 - Auni,k in Eq. 226 is required to 
be zero, i.e., A^m,k is assumed to be constant. Equating Eq. 225 and Eq. 
226 results in the solvable feedforward matrix equation 



r" 

\ 

All ^12 


-A (4) -I) 
i± m 

11 m 

_ H 


A A 

L/I 21 ^22 -1 


H 

m 

K - 

m 


A solution for the feedforward matrices in Eq. 227 exists if the plant has 

no transmission zeroes at 1.0 and no transmission zero equal to an eigenvalue 

of 0 . A numerical solution to Eq. 227 is given in Ref. 5 (the subroutine 
m 

CGTPIF) . The solution is unique if the number of plant controls and plant 
outputs are equal. 

Assuming a snlut±'an with unknown coefficients then determining the values 
of the unknown coefficients by substituting the solution back into the equa- 
tions, as is done for the feedforward control solution, is a technique also 
used in differential equations to determine the homogeneous and particular 
solutions. If the forcing fianction to a differential equation can be repre- 
sented as a model, Eq. 223, with no control input, than the feedforward control 
solution in Eq. 229 can be shown to be the particular solution for that forcing 
function. Assuming the proper form for the solution has been simple for only 
the single-rate discrete- time case. The feedforward control solution for the 
multi-rate case has yet to be determined. 

The feedforward control matrices can be used to relate incremental changes 
in nonlinear trim to incremental changes in the command model similiar to 
results in Eqs. 215 to 217. For a constant command model input, the relationship 
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xs 


'^ll^^o,k“-^o,k-l^ ^2^^o.k"^o,k-l^ 

^,k"^,k-l~ '^21^^o,k~^o,k-l^ ''' ^22^^o,k~^o,k-l^ 


(228) 


(229) 


Small changes in ^ serves the same purpose as small changes in ^ did in 
Eq. 215 and the PIF derivation in Ref. 3, 


u 


u , = u + Au , 


(230) 


X 


A , = X , + Ax , 


(23iy 


The PIF control law is designed in this report to accommodate small changes 


in Au„. 
“-in 
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APPENDIX E 


CONTROL DATA HOLDS ^ COMPUTATION DELAY 
AND TRIM ACCOMMODATION 


This appendix shows how different control data holds and modeled control 
law computation delay affect discrete-time plant model representations. Trim 
accommodation procedures are reviewed. The concepts in this appendix aid in 
understanding the PIF derivation in Chapter IV and command model construction in 
Appendix C. 

A. CONTROL DATA HOLD 

A control data hold is a procedure for constructing in continuous-time 
the control that commands the actuator surface. A zero-order hold, for example 
has the following continuous-time representation 

A^(t) = for ^ ^ ^ (232) 

A triangular data hold has the following representation 

Au(t) = A;^ + (Au^_^^-A^)/At for > t > (233) 


Similiar expressions can be determined for the first-order hold, slewer data 
hold and others. If the intermediate variable 


Av, = 'T-r (Au- , , - Au, ) 
— k At — k+1 — k' 


is defined, then the triangular data hold becomes 


(234) 


Au(t) 

= + (t-tj^)Av^ 

(235) 

AHk+i 

“ AUj^ + At Av^ 

(236) 


The discrete plant representation at the sample points for the triangular 
data hold is obtained by substituting Eq. 235 into Eq. 218, 
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^ Ax(t ) dxGAu + dTGAv, (237) 

k-tx — k o — k o — k 

The equation can be rewritten as 


*^ k+l 

_ 


r ~ 

1 

+ 

■ 

■ * ak + i _ 


_0 


LAUkJ 


— \ 
M 
4-1 
<1 




where 


r , = 


rAt-o Fs , ^ 

/ e ds G 


■p ^ /*C^ ^ ® J r> 

i o ~ e ds G 
X o 


[238] 


(239) 


Equation 238 occurs exactly as shown as a partition when the sampled data 
regulator PIF discrete plant model is formed. The result is altered to 
obtain the desired autopilot zero-order hold mechanization. 

B. COMPUTATION DELAY 


Digital control laws do not necessarily output the control command to 
the actuators at the sample points. There is a delay caused by computation 
and other factors that cause the control command to be released to actuator 
channels between sample points. The control command with computation delay 
can be represented in continuous-time as 

t, ^ > t ^ t, +a 

K+1 k (240) 

t, -Kj > t ^ t, 
k k 


Au(t) = 


Au- 
— k 

^-k-1 


where a is the computational delay. Substituting Eq. 240 into 218 and 
simplifying produces 


^2k+l = + L ^2 


(241) 


where 


r = (f) G 

1 o 




Fs 


dsdT G 


(242) 
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A discrete contro law can accommodate the computation delay, (that is the 
design evaluation and simulation exactly corresponds to the implementation 
response) , by using Eg. 241 as the discrete model representation of the plant. 
Rather than guess at O to perform a design, the PIF autopilot construction 
simply assumes the control is released to the actuators delayed one full cycle. 
Equation 241 with a full cycle delay can be represented as follows 


■^2k+l' 



r" 


+ 

■ 0 ■ 


[243] 

Au, 

1- -I 


_0 

I. 



_AtI 




where Au j^ is the control applied at time Equation 243 (v^ith the u^^ index 

moved up one sample, the implementation is identical) occurs when the discrete 
PIF plant model is simplified. 

C. TRIM ACCOMMODATION 


In typical autopilot operation, the autopilot feeds back an error signal 
to the actuators to drive the error signal to zero. The error signal is form- 
ed by subtracting the sensed signal from a command or trim value. Four popular 
ways to generate trim information in flight are 1) prestore, trim values with 
flight condition, 2) load the measured aircraft states at control initiation 
as the trim values, 3) estimate the trim value of a signal by low-pass filter- 
ing the signal (the error signal is constructed using a "wash-out" filter) 
and 4) estimate trim using integral control. The PIF control law uses the 
second, third and fourth mechanizations. The second mechanization is part of 
the control and command model initialization procedure. The third mechaniza- 
tion is part of the integral control implementation where the rudder command 
error is integrated and fed back. The following simple example demonstrates 
the effect of control error integral feedback. 

Consider a scalar plant 


X = ax + bu 


C244) 


A control law using integral control is formed which feeds back the integral 
of control error from some command value 


u = u^ + X + / (u - u^) 


(245). 


An iraplementable form for Eq. 245 is obtained by taking the derivative and 
simplifying. If the variable 
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z = u - u 

c 


(246) 


is substituted into Eq. 245, the Laplace transfo 2 mi of the control law becomes 


u = z + u 

c 


(247) 



(248) 


Equation 248 is recognized as a wash-out filter. Digitally integrating the 
rudder control command error in a manner similiar to Eq. 245 is used in most 
of the lateral-directional PTF autopilots. The effect is similiar to a yaw 
damper except that all the lateral states fed back are washed out simulta- 
neously- The value of u^ in Eq. 247 for the PIF autopilots is the control 
position at engage. 
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LIST OF SYMBOLS 


In general, matrices are represented by capital letters and vectors are 
underscored; exceptions to these rules are only made when they are contra- 
dicted by standard aerodynamic notation. 


VARIABLE 

A 


ASK 

a 


a 


n 


B 

b 


DESCRIPTION 


Discrete time feedforward matrix 
Fundamental matrix (.continuous-time system) 

Wing reference area 

Gain in ALT SEL command model control system 
Acceleration 

Roll to rudder crossfeed gain 

Component of the earth-relative acceleration in the 
aircraft x-z plane normal to the velocity vector 

Component of the earth-relative acceleration normal to 
the aircraft x-z plane 

Control effect matrix (continuous- time system) 

Wing span 
Bias estimate 


C Control law feedback gains 

C 12 Partial derivative of the nondimens ional coefficient of 

force or moment 1 with respect to the nondimensional 
variable 2 (scalar) 

c Mean aerodynamic chord 

D Control observation matrix 

Drag 

d Perpendicular position from the glideslope to the air- 

craft 


Aerodynamic contact force vector 
_f Vector-valued nonlinear function 

GSK Gain in APR GS command model control system 
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VARIABLE 

g 


H 


HSK 


H 


2 

1 


h 


I 


i 

J 

j 

K 

k 

LK 

LOG 

L 

1 

£ 

M 


m 


N 


n 


DESCRIPTION 


Magnitude of gravitational acceleration vector 
Scalar gain value 

Command observation matrix 

Gain in HDG SEL command model control system 

Euler angle transformation form 
Frame 1 axes to Frame 2 axes 

Vector-valued nonlinear observation function 

Barometric altimeter output 

Identity matrix 
Moment of inertia 

Index integer 

Cost functional matrix 

Gain value 
Index integer 

Gain in APR LOG command model control system 
Gain in APR LOC command model control system 
Aerodynamic moment about the x-axis (scalar) 

Scalar gain value 
Number of commands 

Aerodynamic moment about the y-axis (scalar) 

Cross weighting matrix between states and controls 

Mass of the vehicle 
Number of controls 
Meters 

Aerodynamic moment about the z-axis (scalar) 
Newtons (kg-m/sec”2) 

Number of states 


77 


VARIABLE 

P 

P 

Q 

QZ 

q 

q 

R 

RD 

RLOC 

RMEAS 


r 

S 

s 

T 

t 

UC 

u 

u 

V 

V 


w 

X 

X 


DESCRIPTION 

Riccati matrix in the optimal regulator problem 

Rotational rate about the body x-axis 

State weighting matrix 

Integrator state weighting matrix 

Rotational rate about the body y-axis 

, 1 2 , 

Free stream dynamxc pressure {= 

Control weighting matrix 
Control rate weighting matrix 

Range position from localizer transmitter to touchdown 
point 

Estimate value of range using vertical position and 
glideslope angle 

Rotational rate about the body z-axis 

Feedforward matrix for sensor measurements 

Laplace transform variable 

Thrust 

Time 

Command model control variable 

Body x-axis velocity component 

Control vector 

Velocity magnitude 

Body y-axis velocity component 
Control difference 

Body z~axis velocity component 

Aerodynamic force along the x-axis (scalar) 

Position along the x-axis 


VARIABLE 


DESCRIPTION 


X State vector 

Y Aerodynamic force along the y-axis (scalar) 

y position along the y-axis 

Z Aerodynamic force along the z-axis (scalar) 

z Position along the z-axis 


VARIABLE 

(GREEK) 

Ct 

3 

r 

Y 
6 

a 

^e 

«r 

C 

6 

X 

V 


P 

a 

cr 


DESCRIPTION 

Wind-body pitch Euler Angle (angle of attack) 

Negative of wind-body yaw Euler angle (sideslip angle) 

Discrete time control effect matrix 
Glideslope deviation angle 

Inertial-velocity axis pitch Euler angle (flight-path) 

Aileron deflection 

Elevator deflection 

Rudder deflection 

Damping ratio 

Inertial-body pitch Euler angle 
Eigenvalue 

Euler angle position vector 

Localizer deviation angle 

Eigenvector 
Integrator state 

Air density 

Summation 

Real part of an eigenvalue in radians/sec 
Time delay 
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VARIABLE 

(GREEK) 

T 

$ 

< 1 ^ 

CO 



SUBSCRIPTS 

a 

B 

b 

c 

E 

e 

GS 

H 

i 

j 

k 

1 

m 


n 

o 


DESCRIPTION 

Time constant 

Discrete- time system matrix 

Inertial-body axis roll Euler angle 

Inertial-body axis yaw Euler angle 

Matrix in quad partition matrix inverse 

Frequency in radians/sec 
Imaginary part of an eigenvalue 

Body angular rate vector 

DESCRIPTION 

Aileron 
Body axis 
Bias 

Command value 
Earth- relative axis 
Error quantity 
Glides lope 
Horizontal 

Element index for vectors and matrices 

Element index for vectors and matrices 

Sampling instant index 

Aerodynamic moment about the x-axis 

Aerodynamic moment about the y-axis 
model variable 

Aerodynamic moment about the z-axis 
Nominal value 



SUBSCRIPTS 

DESCRIPTION 

S 

Static pressure 

Accelerometer offset from the eg 

w 

Wind 

X 

Horizontal perpendicular to y and z 

y 

Horizontal perpendicular to x and z 

z 

Vertical perpendicular to x and y 

X 

Aerodynamic force along the x-axis 

y 

Aerodynamic force along the y-axis 

z 

Aerodynamic force along the z-axis 

SUPERSCRIPTS 

DESCRIPTION 

E 

Earth (inertial) axis 

T 

Transpose of matrix 

W 

Wind axis 

-1 

Inverse of matrix 

* 

Steady state 
Trim 

PUNCTUATION 

DESCRIPTION 

(*) 

Derivative of quantity with respect to time 

(J 

Vector quantity 

3( )/3( ) 

Partial derivative of one variable with respect to 
another 

A( ) 

Perturbation variable 

( )* 

Star trajectory 

(") 

Estimated quantity 

Discrete cost fxanction weighting matrix 

00 

Infinity 
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PUNCTUATION 


DESCRIPTION 


f Integral 

[• ] Important Equation 

(• ) Equation 

(~) Matrix equivalent to vector cross product. 

Specifically, if x is the three-dimensional vector 



and the cross product of x and ^ is equal to the 
product of the matrix x and the vector f, 

X X f = xf 


ACRONYM 

CORRESPONDING PHRASE 

ARA 

Avionics Research Aircraft 

eg 

Center of gravity 

DME 

Distance Measuring Equipment 

GA 

General Aviation 


GATOR 

General Aviation Terminal Area 

Operation Research 

ILS 

Instrument Landing System 


kt 

Knot 


MAX 

Maximvim 


MLS 

Microwave Landing System 


NASA 

National Aeronautics and Space 

Administration 

PI 

Proportional- Integral 


PIF 

Proportional- Integral- Filter 


rad 

Radians 


sec 

Seconds 
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ACRONYM 


VALT 


CORRESPONDING PHRASE 


VTOL Approach and Landing Technology 
Vertical Takeoff and Landing 


VTOL 



TABLE 1 NAVION ACTUATOR CHARACTERISTICS 


CONTROL 

DISPLACEMENT 
LIMIT, deg 

RATE LIMIT, 
deg/sec 

TIME CONSTANT, 
T, (sec) 

MAXIMUM SPECIFIC 
MOMENT (IAS=53.8m/s 
(105kt)) 






AILERON 

±20.0 

70.0 

0.17 

9.2 rad/sec^ 





roll acceleration 

ELEVATOR 

-30.0 

70.0 

0.17 

9.9 rad/sec^ 


+20.0 



pitch acceleration 

RUDDER 

±20.0 

70.0 

0.17 

4.2 rad/sec^ 


TABLE 2 SENSOR NOISE AND BIAS CHARACTERISTICS 
FOR THE NAVION RESEARCH AIRCRAFT 


SENSOR 

NOISE STANDARD DEVIATION 

BIAS STANDARD DEVIATION 




Lateral 

Accelerometer 

0.3 m/s^ 

0. 6 m/ s^ 

Vertical 

Accelerometer 

2 

0.6 m/s 

0. 6 m/ 

Rate Gyro 

0.2 deg/sec 

0.2 deg/sec 

Attitude Gyro 

0.1 deg 

0. 2 deg 

Heading Gyro 

0.1 deg 

0.5 deg 

Indicated Airspeed 

0.05 m/s (0.15 fps) 

0.15 m/s (0.5 fps) 

Barometric Altimeter 

3.0 m (10 ft) 

15.0 m (50 ft) 

Localizer /Glides lope 
Angular Deviation 

0.02 deg 

0.012 deg 

Outer Marker 
Indication 

— 

500.0 m (1640 ft) 
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TABLE 3 ROLM 1666 CHARACTERISTICS 


DESCRIPTION 

MEMORY 


General purpose 16 - bit minicomputer 
designed to MIL-E-5400 specifications 

65,536 words of 1 ysec ferrite core 


EXECUTION TIME - (Time in ysec, register to register 

operations) 


INSTRUCTION 


FIX ED POIN T FLOATING POINT 

■ (16 bit) ” (32 bit) 


Add 


Multiply 

5.2 

Divide 

9.2 

Load, Store 



.0 1.8 - 4.8 

- 5.4 3.6 - 4.8 

- 9.6 8.0 - 8.8 

.0 4.8 



TABLE 4 SAMPLE SUMMARY OF COMMERCIALLY AVAILABLE AUTOPILOTS 



lOOOA 

THIS REPORT 
PIF 

















TABLE 5 STATES AND CONTROLS FOR AUTOPILOT MODE DESIGNS 


BETA HOLD 




MODEL 


i 

u 

-TH 

= 

m 



T 

X 

— m 

= 

in 



T 

= 

m 

(t> ] 


T 

X 

= 

[6 

r p (j)] 


T 

_u 

= 


6 ] 
r 


T 

Z 

= 

[3 

<1>] 

ROLL SEL 




MODEL 


T 

u 

— m 

= 


5 ] 

rm 


T 

X 

— m 

= 


■s ] 

rm 


T 

= 

m 

(a, (J) + <5 )] 

k m rm 


T 

5C 

= 

[a.y 

r p (f)] 


T 

u 

= 

'«a 



T 

Z 

= 



HDG SEL 




NONLINEAR MODEL 


T 

u 

—m 


[L 

m 

6 ] 
rm 


T 

X 

— m 

- 

m 

(p ] 

m 


T 

= 

m K m rm 


T 

X 

= 

[ay 

r p 4) i|;] 


T 

_u 

= 


<S ] 

r 


T 

Z 

= 

[4> 

6 ] 
r 
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HDG SEL 


LINEAR MODEL 



A T 

Au = 
— m 

[A(|) A6 ] 

m rm 


T 

Ax^ = 

wj 


II 

-^k + “rm' 


II 

[ Aa ^ Ar Ap A(j) ] 


A T 

Au = 

[A6 A6 ] 

a r 


II 

[A(})-Ai|; -a^ Aip + 6^] 

ALT SEL 


MODEL 


u = 
— m 

m 


T 

X = 

—m 

[z Z ] 

m m 



m 


II 

[V a q 0 z] 

n 


u = 

[y 


y = 

[z] 



LINEAR MODEL 


Au = 
— m 

m 


Ax = 
“in 

m 


% = 

m 


II 

[AV Aa Aq A0 Az] 

n 


Au = 

[A6 ] 
e 


Ay = 

[Az] 




TABLE 5 (CONTINUED) STATES AND CONTROLS FOR AUTOPILOT MODE DESIGNS 


PITCH 


ARP GS 


SEL 


MODEL 


u = 


m 


m 


X = 


m 


m 




m 


in 


X = [V 


q 6] 


u = 


[«el 


y - [ej 


NONLINEAR MODEL 


m 


Ttl 


X = 


m 


'"ml 


y„ - 


m 


m 


X = 


[V a q 6 d] 
n 


u = 


[«el 


[d] 


LINEAR MODEL 


Au = 

[AX ] 

m 

m 

Ax = 

[Ad ] 

la 

m 


[Ad ] 

m 

m 

Ax^ = 

[AV 


ix^ = [AV Aa Aq A9 Ad] 

— n * 



TABLE 5 (CONTINUED) STATES 


APR LOG 


APR LOCI 


APR LOCR 


APR LOCP 


APR LOG 


APR LOCI 


APR LOCR 


APR LOCP 




TABLE 6 AUTOPILOT MODEL DESIGN PARAMETERS 


AUTOPILOT 
HDG SEL 

ALT SEL 
APR GS 


APR LOG 


APR LOCR 
APR LOCP 


EQUATION 

PARAMETERS 

VALUE 

UNITS 





150 

HSK 

1.0 



150 

HSK 2 

2.0 

— 

155 

^maxl 
Tclose 
, maxi 

0.262 

radians 

— 

0.157 

radians 

153 

0.0873 

radians 

145 

T 

10.0 

— 

154 

914.4 

m 

162 

ASK. 

ASK^ 

..close 

0.8 

— 

162 

2.0 

— 

— 

6.0 

m 2 

165 

fmaxl 

0.2286 

m/ s 

167 


2.53 

m/s 

166 

a 

max 

0.24 

radians 

177 

OUTMARK 

9879.0 

m 

— 

RMIN 

100.0 

m 



r- 

0.0436 

radians 

182 

gmax 

-0.0004 

ft"^ 

182 

82 

0.07 

— 

— 

S 3 

0.07 

— 

— 

S 4 

0.04 

— 

184 


0.0523 

radians 

— 

0.0 

radians 

^ 

^lUdLX 

-0.08726 

radians 



0.000524 

radians 

— 

max 

-0.00037 

m“^ 

— 


0.112776 

— 

— 

RKMAX 

0.08 

— 

200 

RLOC 

2205.0 

m 

— 

77 

0.0122 

radians 

201 

.max 

X.. 

- 0.0002 

ft"l 

201 

ti 

0.2 

— 

— 


0.2 

— 

— 

^4 

0.11 


203 

8 

9.81 

/ ^ 
Ttl/S 

203 

LOG 2 

1.2 

— 

— 

’^cmax 

0.523 

radians 

— 

tmax 

0.175 

radians 


*^max 

0.00873 

radians 

191 

k, 

y 

3.0 



191 

0.017 

— 

192 

k 

y 

0.04 

— 


AUTOPILOT MODE 


* 

STATE 


OR 

COOTROL 

ALT SEL 
(DESIGN) 

ALT SEL 

(LOWER 0 WEIGHT) 

PITCH SEL 

APR GS 






A6e 

0.0 

0.0 

0.0 

0.0 


— 

— 

— 

— 

A6^ 

— 

— 

— 

— 

Au^" 

0.0 

0.0 

0.0 

0.0 

Aw 

0.0 

0.0 

0.0 

0.0 

Aq 

0.0 

0.0 

0.0 

0.0 

A0 

11.0 

5.0 

8.0 

12.5 

Az 

0.5 

0.5 

0.0 

0.0 

Ad 

— 

— 

— 

0.5 

Av 

— 

— 

— 

— 

Ar 

— 

— 

— 

— 

Ap 

— 

— 

— 

— 

A4> 

— 

— 

— 

— 

A^ 


— 

— 

— 

Ay 


— 

— 

— 

^Clong 

0.25 

0.25 

3.0 

0.25 

^^latl 

— 

— 


— 

^|lat2 


— 

— 

— 

A6e 

7.0 

7.0 

7.0 

7.0 

A^a 

— 

— 

— 

— 

A6r 

— 


— 

— 

(Az or Ad) 

1.0 

1.0 

0.0 

1.0 













TABLE 7 (CONTINUED) SQUARE ROOT OF Q AND R DIAGONAL 
WEIGHTS IN THE PIF COST FUNCTIONS 


STATE* 

OR 

CONTROL 

AUTOPILOT MODE 

APR LOCI 

APR LOCR 

APR LOCP 




A6 

0.1 

0.1 

0.1 

A6^ 

0.1 

0.1 

0.1 

Au^ 

Aw 

Aq 

A6 

Az 

Ad 

Av 

___ 

,-r-i 


0.0 

0.0 

0.0 

Ar 

11.0 

11.0 

11.0 

Ap 

0.0 

0.0 

0.0 

Acj) 

10.0 

10.0 

10.0 


14.0 

11.0 

14.0 

Ay 

0.0 

0.0 

0.0 

/Ay 

2.5 

1.0 

1.0 


0.15 

2.2 

2.2 


7.0 

7.0 

7.0 

A6^ 

r 

7.0 

7.0 

7.0 


* units are degrees and meters 
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TABLE 8 OPEN- AND CLOSED-LOOP s-DOMAIN EIGENVALUES 


OPEN LOOP 


DYNAMIC 

MODE 


DYNAMIC 

MODE 


SHORT PERIOD 
PHUGOID 




ALT SEL (DESIGN) 

ALT SEL (LOWER 0 WEIGHT) 

T, 

rad/sec - sec 

C. T, 

rad/sec - sec 


3.0 0.78 

0.25 0.025 


PITCH SEL 


3.07 0.69 


1.58 0.99 

0.25 0.77 


3.07 0.77 


1.12 0.81 

0.36 0.80 



3.04 

0.74 

^ 

1.35 

0.91 

— 

— 

— 

8.3 

— 

— 

2.7 

OPEN LOOP 


. 

T, 

rad/sec 

sec 

1 — _ __■ 

2.07 

0.22 

— 

— 

— 

0.16 



-31.0 

-.^OO 

APR LOCI 


3.21 0.67 


0.24 0.77 


/0.74t 

y0.50/ 



DYNAMIC 

MODE 


DUTCH ROLL 
ROLL 
SPIRAL 
HEADING 
DUTCH roll-6 
16 -V ' 


DYNAMIC 

MODE 


ROLL HOLD 




HDG SEL 


5 , 


rad/sec 



APR LOCR 



APR LOCP 



T, 

rad/ sec - 

sec 


DUTCH ROLL-d 

2.9 

0.47 

— 

2.84 

0.48 

— 

2.9 

0.47 

— 

ROLL ' 

— 

— 

0.16 

— 

— 

0.16 

— 

— 

0.16 

/y-v 

0.29 

0.58 

— 

1 0.36 

0.56 

— 

0.35 

0.61 

— 

6a-<' 

1.73 

0.69 

— 

1.72 

0.69 

— 

1.77 

0.69 

— 

y 




1.05 

2.00 

(0-.60 

0.99)^ 



— 

1.07 

3.08 

M 

— 

— 

5.9 

— - 

— 

5.6 

— 

— 

6.04 


The mode split into two real roots. 

The roots combined to form a complex pair. 
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TABLE 9 ELEVATOR CONTROL GAINS FOR THE LONGITUDINAL AUTOPILOTS 




AUTOPILOT 

MODES 


STATE* 

ALT SEL 

ALT SEL 

PITCH SEL 

APR GS 


(DESIGN) 

(LOWER 0 WEIGHT) 




1 

Ve 

0.719 

0.79 

0.753 

0.704 

V 





r 

V 

-0.025 

-0.018 

-0.00928 

-0.027 

1 

a 

-0.021 

-0.015 

-0.011 

-0.022 

q 

-0.428 

-0.28 

-0.294 

-0.466 

0 

-2.2 

-1.48 

-1.42 

-2.34 

z 

0.0139 

0.012 

0.0 

0.0 

r 

^ 

__ 





P 






— 

d 

— 

— 

— 

-0.0144 

y 

flong 

^latl 

Hat2 

mlong 

0.000176 

0.000183 

-0.038 

-0.000175 

-0.0491 

-0.029 

1.08 

2.49 

'^mlatl 

^mlat2 

— 

— 

— 

— 


* units 


are radians and meters 




TABLE 10 CONTROL GAINS FOR THE ROLL SEL AND HDG SEL AUTOPILOTS 



AUTOPILOT MODES 

STATE* 


AILERON 

RUDDER 1 


ROLL SEL 

HDG SEL 

ROLL SEL 

HDG SEL j 



1 

V 

— 

— 

— 



e 

V 

0.65 

0.69 

0.021 

0.061 

a 

V 

r 

0.020 

0.063 

0.76 

0.715 

^T 





z 

q 

0 

— 

— 

— 

— 

z 

a 

-0.053 

0.012 

0.178 

0.537 

r^ 

-0.126 

0.25 

-0.74 

-2.0 

P 

-0.31 

-0.26 

0.011 

0.087 

<P 

-2.16 

-1.51 

-0,53 

-0.44 

4> 

d 

— 

0.50 

— 

-5.9 

y 

^long 

^latl 

^lat2 

-0,068 

-0.014 

-0.022 

0.032 

0.019 

0.028 

-0,512 

0.013 

mlong 

%latl 

^mlat2 

1.5 

1.15 

0.67 

1.28 

-0.35 

- ■ 

-1.19 

0,81 

-0.59 


* units are radians and meters 









AUTOPILOT MODE 


STATE* 

APR LOCI 

APR LOCR 

APR LOCP 





V 

— 

— 

— 

e 

V 

0.69 

0.69 

0.69 

a 

V 

0.032 

0.048 

0.074 ! 





V 



— — — 

T 




a 

— — — 

“ 

“““ 

z 




q 



“““ 

0 

— 

— 

— 

z 

— 

— 

— 

a 

0.18 

0.278 

0.39 


-0.26 

-0.356 

-0.38 

P 

-0.25 

-0.245 

-0.25 


-1.60 

-1.55 

-1.54 


-2.15 

-2.95 

-4.19 

d 

— 

— 

— 

y 

-0.019 

-0.014 

-0.026 

E 

?long 

;iati 

^lat2 

-0.021 

0.000798 

0.0075 

-0.00026 

0.0269 

0.0012 

long 

""latl 

1.57 

1.74 

-1.25 

U , r, 

0.0 

-1.53 

2.05 

mlat2 









TABLE 12 RUDDER CONTROL GAINS FOR THE APPROACH AUTOPILOTS 




AUTOPILOT MODE 


STATE* 

APR LOCI 

APR LOCR 

APR LOCP 

1 




V 




e 




V 

a 

0.000277 

-0.000674 

-0.00068 

V 

r 

0.827 

0.821 

0.824 





T 




a 






z 




q 

— 

— 

— 

0 

— 

— 

— 

z 

— 

— 

— 

a 

y 

0.73 

0.706 

0.67 

r 

-1.22 

-1.24 

-1.18 

P 

-0.027 

-0.021 

-0.027 

(t> 

-0.61 

-0.66 

-0.65 


-5.54 

-5.54 

-5.45 

d 


— 

— 

y 

-0.0224 

-0.031 

-0.028 

y 

long 

— 

— 

— 

^latl 

0.0227 

0.0123 

0.00967 

^lat2 

-0.000237 

-0.00171 

-0.000947 

^mlong 

— 

— 

— 

^mlatl 

1.41 

1.59 

1.16 

%lat2 

0.0 

0.310 

1.72 


* units are radians and meters 





TABLE 13 MATRIX ELEMENTS FOR THE PIF AUTOPILOT MODES 


ALT SEL 


PITCH SEL 


APR GS 


ROLL SEL 


■ 0 " 


" 0 

0 


0 

0 


0 

0 


1.0 

1 . 0 _ 

1 


‘ 0 “ 


■ 0 o' 

0 


0 0 

0 


0 0 

0 


1.0 0 

l-0_ 



HDG SEL 


" 0 ■ 
0 
0 
0 

1.0 


APR 

LOCI 

APR 

LOCR 

APR 

LOCP 

0 

0 ■ 


- 0 

0 “ 


■ 0 

0 ■ 

0 

0 


0 

0 


0 

0 

0 

0 


0 

0 


0 

0 

0 

0 


0 

0 


0 

0 

0 

1.0 


0 

1.0 


0 

1.0 

1.0 

0 . 


.1.0 

0.31. 


.1.0 

0 _ 


TABLE 14 NAVION AIRPLANE PARAMETERS 


Wing area, A 

Mean Aerodynamic Chord, c 
Wing Span, b 
Gross Mass, WT 

ix 

-z 

^xz 

17.112 m^, (184 ft-) 

1.74 m (5.7 ft) 

10.17 m (33.38 ft) 

- 1540.6 kg (3400 lb) 

- 1742.33 kg-m2 (1284.08 slug-ft^) 

- 3762.4 kg-m2 (2772.86 slug-ft2) 

- 4389.1 kg-m2 (3234.72 slug-ft2) 

0.0 

Control Surface 

Deflection, deg 

Flaps 

, 0 to -40 

Elevator 

- -19 to 29 

Aileron 

- -18 to 19 

Rudder 

- -20 to 25 




TABLE 15 FLIGHT CONDITIONS FOR AERODYNAMIC DATA 



Ref. 31, 
Condition I 

Ref. 31, 
Condition II 

Ref. 38, 

Altitude 

1524.0 m 
(5000 ft) 

1524.0 m 
(5000 ft) 

1980.0 m 
(6500 ft) 

Velocity 
(true airspeed) 

73.2 m/sec 
(2400 ft/sec) 

44.0 m/sec 
(144.0 ft/sec) 

54.0 m/sec 
(176.0 ft/sec) 

Flaps 

0 deg 

20 deg 

10 deg 
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TABLE 16 LONGITUDINAL AERODYNAMIC PARAMETERS 


Parameter 

Condition I, 

Condition II, 



Ref. 31 

Ref. 31 

Ref. 38 



Design Values 

Analog Match 


0.262 

-1.37 

— 

^Xa 

-A. 33 

-4.86 

-4.6 

^Zq 

-15.9 

-27.13 

0.0 

^6e 

0.511 

0.52 

0.29 


-0.77 

-0.84 

0.0 


-6.5 

-6.0 

-0.64 

% 

•18.1 

-16.4 

-14.4 


1.42 

1.55 

0.96 


0.0055 rad 

0.105 rad 

0.036 rad 

‘lo 

0.0 rad 

0.0 rad 

0.0 rad 


0.4 m/sec(1.3 ft/sec) 

4.6 m/sec(15.2 ft/sec) 

1.1 m/sec (3.7 ft/sec) 


73 m/sec(240.0 ft/sec) 

44 m/sec(144.1 ft/sec) 

55 m/sec (185 ft/sec) 

“o 

0.0055 rad 

0.105 rad 

0.02 rad 

So 

0,0015 

0.0015 

— 

So 

-0.27 

-0.75 

— 

^Md 

0.0 

0.0 

— 

SSj 

0.0061 

0.0061 

— 

g 

9.8 m/sec^ 

9.8 m/sec^ 

9.8 m/sec^ 

p 

.0.00912 kg/m^ 
(0.00205 slugs/ft^) 

0.00912 kg/m^ 
(0.00205 slugs/ft ) 

0.00912 kg/m^ 
(0.00205 slugs/ff^) 



TABLE 17 LATERAL AERODYNAMIC PARAMETERS 


Parameter 

Condition I, 
Ref. 31 

Condition II, 
Ref. 31 
Design Values 

Ref. 38 
Analog Match 

S 

0.0 

0.0 

0.0 

Sr 

0.0 

0.0 

0.0 

Sb 

-0.6 

-0.74 

-0.92 

"?6r 

-0.33 

-0.143 

-0.21 


-0.026 

-0.023 

-0.025 


-0.07 

-0.053 

-0.086 

% 

-0.49 

-0.53 

-0.52 

^Lr 

0.11 

0.114 

0.093 

^L6a 

0.154 

0.16 

0.16 


0.0 

0.0 

-0. 025 

Sb 

0.073 

0.080 

0.103 

Sp 

-0.04 

-0.147 

-0.049 

Sir 

-0.09 

-0.12 

-0.11 


-0.004 

-0.0015 

-0.0039 

<i.a. 

0.063 

0.075 

0.11 

V 

o 

0.0 m/sec 

0.0 m/sec 

0.0 m/sec 


0.0 rad 

0.0 rad 

0.0 rad 

Po 

0.0* rad 

0.0 rad 

0.0 rad 


0.0 rad 

0.0 rad 

0.0 rad 

'i>o 

0.0 rad 

0.0 rad 

0.0 rad 

Sfo 

0.0 

0.0 

0.0 

^Lo 

0.0 

0.0 

0.0 

So 

0.0 

0.0 

0.0 





WALLOPS FLIGHT CEfs^TER FACILITIES 



FIGURE 1 WALLOPS FLIGHT TEST FACILITIES 
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FIGURE 2 BASIC PIF CONTROL LAW 
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FIGURE 3 COMPARISON OF OPEN-LOOP LINEAR MODEL RESPONSE 
WITH AIRCRAFT RESPONSE FOR A COMMAND IN RUDDER 



FIGURE 3 (CONTINUED) Comparison of Open-Loop Linear Model 

Response with Aircraft Response for a Command in Rudder 



TIME ( SEC) 



FIGURE 4 COMPARISON OF OPEN-LOOP LINEAR MODEL RESPONSE 

WITH AIRCRAFT RESPONSE FOR A COMMAND -IN ELEVATOR 
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ARA CONFIGURATION 



AIRCRAFT SYSTEM 


STANDARD 

PANEL 

DISPLAYS 


VISUAL & 

INERTIAL 

CUES 


n 


STANDARD 

COCKPIT 

CONTROLS 


CONTROL 
. SURFACES & 
THROHLE 


RESEARCH 

AIRCRAFT 


CONTROL SYSTEMS 




DARES YSTEMI 

r 

SYSTEM 


MOTION 

SENSOR 


AIRCRAFT 


I MOTION 


J 


EXPERIMENTAL 
COCKPIT 
CONTROLS 


EXPEF IMENTAL 
DISPIAYS 


VISUAL & 
INERTIAL 
CUES 
I 


FIGURE 5 ARA CONFIGURATION 


AVIONICS RESEARCH AIRCRAFT CONTROL/DISPLAY PANEL 


EVALUATION PILOT 
STATION 


SAFETY PILOT 
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ANNUNCIATORS [ 

© 
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-®c 
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) 

D o 
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V 

/I 

i 
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Ui 

J 


MODE SWITCHES 
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( 

[VARIABLE 

1 

CIRCUIT 
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1 BREAKERS 

IPGTENTIOMETERS 


L, J 

L > 


ENGINE INSTRUMENTS 


GEAR VSS PANEL 


• ^'CRT: 


im 

±iii 


KEYROARD/j j| . 
^DISPLAY_n|^ 


SWITCH I VARIABLE -STABILITY 

MAIN SYSTEM (VSS) 

POWER 1 POTENTIOMETERS 


FIGURE 6 AVIONICS RESEARCH AIRCRAFT CONTROL /DISPLAY PANEL 
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FIGURE 7 DETAILED PIF CONTROL LAW BLOCK DIAGRAM 















MAIN EXECUTIVE 



FIGURE 8 FLOW DIAGRAM FOR THE MAIN EXECUTIVE 
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FIGURE 9 (CONTIITOED) FLOW DIAGRAM FOR THE SUBROUTINE CONTROL 
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FIGURE 11 FLOW DIAGRAM FOR THE SUBROUTINE FPIFl 
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PIF CGT 
PROGRAM 


FIGURE 12 PIF DESIGN SEQUENCE 
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FIGURE 14 SPIRAL DESCENTS USING THE HDG SEL AND ALT SEL PIF 
AUTOPILOTS 
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FIGURE 15 A 45 DEG HEADING CHANGE USING THE HDG SEL AUTOPILOT 
AND AN a PREFILTER WITH A 0.01 SEC TIME CONSTANT 

y 
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FIGURE 16 A A5 DEG HEADING CHANGE USING THE HDG SEL AUTOPILOT 
AND AN a PREFILTER WITH A 0.10 SEC TIME CONSTANT 
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AILERON (DEG) , ROLL ( DEG) BETA ( DEG) 


6 



0 11 22 33 
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FIGURE 16 (CONTINUED) A 45 DEG HEADING CHANGE USING THE HDG SEL 

AUTOPILOT AND AN a PREFILTER WITH A 0.10 SEC 
TIME CONSTANT ^ 
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25 36 ^7 58 63 

TIME (SEC) 


FIGURE 17 A 45 DEG HEADING CHANGE USING THE HDG SEL AUTOPILOT 
AND AN a PREFILTER WITH A 0.22 SEC TIME CONSTANT 



ELEVATOR (DEG) 










FIGURE 19 A 30.5 m (100 FT) ALTITUDE DESCENT USING THE ALT SEL 

AUTOPILOT, NO THROTTLE ADJUSTMENT AND THE DESIGN GAIN SET 
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FIGURE 20 (CONTINUED) ALT SEL LINEAR AIRCRAFT MODEL RESPONSE SIMULATION 
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CONTROL RESPONSE PIRCRBFT RESPONSE ftIRCRfiFT RESPONSE 








R[RCRflFT RESPONSE 
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FIGURE 22 (CONTINUED) ROLL SEL LINEAR AIRCRAFT MODEL RESPONSE SIMULATION 
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FIGURE 23 (CONTINUED) APR LOCI LINEAR AIRCRAFT MODEL RESPONSE SI>1ULATI0N 
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FIGURE 24 (CONTINUED) APR LOCR LINEAR AIRCRAFT MODEL RESPONSE SIMULATION 
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FIGURE 25 (CONTINUED) APR LOCP LINEAR AIRCRAFT MODEL RESPONSE SIMULATION 
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FIGURE 27 A 152.4 m (500 ft) ALITUTDE ASCENT USING THE ALT SEL AUTOPILOT, THROTTLE 
ADJUSTMENTS, LOW PITCH WEIGHT GAIN SET AND HEAVY TURBULENCE CONDITION 
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FIGURE 28 A 152.4 m (500 ft) ALTITUDE ASCENT USING THE ALT SEL AUTOPILOT, THROTTLE 
ADJUSTMENTS, LOW PITCH WEIGHT GAIN SET AND MEDIUM TURBULENCE CONDITIONS 
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FIGURE 29 (CONTINUED) A 152.4 m (50,0 ft) ALTITUDE ASCENT USING THE ALT SEL AUTOPILOT, 

THROTTLE ADJUSTMENTS, DESIGN GAIN SET AND LIGHT TURBULENCE 
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TIME (SEC) 


FIGURE 30 (CONTINUED) A 152.4 m (500 ft) ALTITUDE ASCENT USING THE ALT SEL AUTOPILOT 

NO THROTTLE ADJUSTMENTS, DESIGN GAIN SET AND LIGHT TURBULENCE 
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FIGURE 32 A 152.4 m (500 ft) 1>1ANUAL ALTITUDE DESCENT WITHOUT USING 
LIGHT TURBULENCE 
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FIGURE 34 A 5 deg PITCH CHANGE USING THE PITCH SEL AUTOPILOT 
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FIGURE 36 APR GS GLIDESLOPE CAPTURE AND TRACK WITH OPEN-LOOP RANGE ESTIMATION AND 
THE AIRCRAFT CROSSING THE OUTER MARKER 
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FIGURE 37 APR GS GLIDESLOPE CAPTURE AND TRACK WITH CLOSED-LOOP RANGE ESTIMATION AND 
THE AIRCRAFT CROSSING THE OUTER MARKER 
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FIGURE 38 APR GS GLIDESLOPE CAPTURE AND TRACK WITH CLOSED-LOOP RANGE ESTIMATION 
AND THE AIRCRAFT MISSING THE OUTER MARKER 
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FIGURE 38 (CONTINUED) ARP GS GLIDESLOPE CAPTURE AND TRACK WITH CLOSED-LOOP RANGE 

ESTIMATION AND THE AIRCRAFT MISSING THE OUTER MARKER 
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FIGURE 39 APR LOCR LOCALIZER CAPTURE AND TRACK FOR A 45 deg INTERCEPT ANGLE 
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FIGURE 40 APR LOCP LOCALIZER CAPTURE AND TRACK FOR A 45 deg INTERCEPT ANGLE 
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FIGURE 41 APR LOCI LOCALIZER CAPTURE AND TRACK FOR A 45 deg INTERCEPT ANGLE 
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FIGURE 43 (CONTINUED) 
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FIGURE 45 


GLIDESLOPE GEOMETRY 


X - Axis 



FIGURE 46 LOCALIZER GEOMETRY 
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